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DOMAINS OF ANALYTICITY OF LINDSTEDT EXPANSIONS OF KAM 
TORI IN DISSIPATIVE PERTURBATIONS OF HAMILTONIAN 

SYSTEMS 

RENATO C. CALLEJA, ALESSANDRA CELLETTI, AND RAEAEL DE LA LEAVE 


Abstract. Many problems in Physics are described by dynamical systems that are con¬ 
formally symplectic (e.g., mechanical systems with a friction proportional to the velocity, 
variational problems with a small discount or thermostated systems). Conformally symplec¬ 
tic systems are characterized by the property that they transform a symplectic form into a 
multiple of itself. The limit of small dissipation, which is the object of the present study, is 
particularly interesting. 

We provide all details for maps, but we present also the (somewhat minor) modifications 
needed to obtain a direct proof for the case of differential equations. We consider a family 
of conformally symplectic maps defined on a 2c?-dimensional symplectic manifold A4 
with exact symplectic form fl; we assume that satisfies /* = A(e)D. We assume 

that the family depends on a d-dimensional parameter /r (called drift) and also on a small 
scalar parameter e. Furthermore, we assume that the conformal factor A depends on e, in 
such a way that for e = 0 we have A(0) = 1 (the symplectic case). We also assume that 
A(e) = 1-1- Q;e“ -I- 0(|e|“+^), where a S Z+, a S C\{0}. 

We study the domains of analyticity in e near e = 0 of perturbative expansions (Lindstedt 
series) of the parameterization of the quasi-periodic orbits of frequency to (assumed to be 
Diophantine) and of the parameter Notice that this is a singular perturbation, since any 
friction (no matter how small) reduces the set of quasi-periodic solutions in the system. We 
prove that the Lindstedt series are analytic in a domain in the complex e plane, which is 
obtained by taking from a ball centered at zero a sequence of smaller balls with center along 
smooth lines going through the origin. The radii of the excluded balls decrease faster than 
any power of the distance of the center to the origin. We state also a conjecture on the 
optimality of our results. 

The proof is based on the following procedure. To find a quasi-periodic solution, one solves 
an invariance equation for the embedding of the torus, depending on the parameters of the 
family. Assuming that the frequency of the torus satisfies a Diophantine condition, under 
mild non-degeneracy assumptions, using a Lindstedt procedure we construct an approximate 
solution to all orders of the invariance equation describing the KAM torus; the zeroth order 
Lindstedt series is provided by the solution of the invariance equation of the symplectic case. 
Starting from such approximate solution, we implement an a-posteriori KAM theorem to 
get the true solution of the invariance equation, and we show that the procedure converges. 
This allows also the study of monogenic and Withney extensions. 
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1. Introduction 


Many problems of physical interest are described by models given by Hamiltonian systems 
with a small dissipation. 

Some important particular cases of systems with dissipation are the following: 

• Hamiltonian systems with a dissipation proportional to the velocity describing, e.g., 
problems of Celestial Mechanics - see [CTO] : 

• Euler-Lagrange equations of exponentially discounted systems - very natural in fi¬ 
nance, when inflation is present and one wants to minimize the cost in present money 
- see |Ren88llTSMrillDFTZ14] : 

• Gaussian thermostats - see [DMOGl IWLQSj . 

In the examples above, it was discovered that there is a nice geometric structure, namely 
that the natural symplectic form is transformed into a multiple of itself by the dynamics. 
Systems that have this geometric property are referred to as conformally symplectic. Besides 
their applications to physical problems, conformally symplectic systems were considered on 
their own in differential geometry ( |Bann21 AgrlO| ). 

This geometric structure has important consequences for the dynamics (see Section [TTT]) . 
Notably for our purposes, there is a KAM theory with an a-posteriori format and a systematic 
way of obtaining perturbative expansions, see |CCdlL13^ . 

In this paper, we study the analyticity properties of the parameterization of the quasi- 
periodic solutions and of the drift parameter of conformally symplectic systems. Perturbative 
expansions to all orders are easy to obtain and have been considered for a long time (as a 
matter of fact, we will develop also very efficient methods of computation of perturbative 
expansions). 

Notice that adding a dissipation to a Hamiltonian system is a very singular perturbation. 
We expect that a Hamiltonian system admits quasi-periodic solutions with many frequencies. 
On the other hand, a system with a positive dissipation - even if extremely small - leads to 
the creation of attractors which have few - or even none ! - quasi-periodic solutions. This 
is why one has to consider external parameters such as the drift. 

Hence, we do not expect that the asymptotic expansions converge and therefore the pa¬ 
rameterization functions will not be analytic in e in any ball centered at the origin. A 
common strategy in the literature of singularly perturbed problems has been to device as¬ 
ymptotic expansions and show that, even if they do not converge, they can be resummed. See 
|Bal941 IHar49[ IBGM96] for general mathematical treatments and |LGZJ90l IAFG901 IKSFOl] 
for surveys of applications of resummation techniques to concrete problems in Physics. For 
singular perturbation expansions in dynamical problems similar to the ones in this paper, 
there have been quite a number of studies (see, e.g., |GM961100021 lOBDOO] lOenlOl 10005] 
IGF0131 IGF014| b In this work, we follow a different approach from resummation. By using 
an a-posteriori theorem, we show (see Theorem [T^ that the parameterization and the drift 
are analytic in e in a domain (called Q and defined precisely in (12.bh ) obtained by removing 
from a ball centered at the origin a sequence of (much smaller) balls with centers in smooth 
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curves going through the origin. The radii of these balls decrease very fast (faster than any 
power of the size of the distance to the center of the excluded ball) as the centers of the 
excluded balls go to zero. Hence, even if the domain does not contain a ball centered at the 
origin, it is hard to distinguish it from a full ball (see Figured]). 

We emphasize that we only prove rigorously that ^ is a lower bound for the analyticity 
domain. The results can be improved in many ways (see Section flO.dh . On the other hand, 
in Section [9] we conjecture that the domain Q is essentially optimal in the sense that for 
a generic system none of the excluded balls can be hlled completely (we present heuristic 
arguments for the conjecture and a proof of a weaker result (Proposition |23|) that shows that, 
given any ball, one can hnd systems where the analyticity domain does not cover the ball). 
We note that the domains we hnd have only sectors centered at the origin with very small 
apertures, say tt/o, where a is the leading exponent in the change of the conformal factor. 

The conjecture that there are no sectors of analyticity with aperture bigger than vr/a has 
consequences for resummability methods and asymptotic expansions, since uniqueness of the 
function given by the expansion and Borel summation often requires that there are sectors 
of angular aperture of about tt, see |Sok80] (the well known Cauchy functions exp( —1/e^) 
provide examples of non-trivial functions with trivial asymptotic expansions in domains of 
angular aperture less than tt). 

We also recall that, since the parameterization is a function, the notion of analytic func¬ 
tions taking values on parameterizations requires that we specify a Banach space in which 
the parameterization lies. This is done in Section [2j We anticipate that these spaces for the 
parameterization are spaces of analytic functions from a complex extension of the torus to 
the phase space. 

The method of proof of this paper is to formulate a functional equation for the parame¬ 
terization of the torus and the drift expressing that the torus is invariant for the map with 
the adjusted parameters. First, we show (part A of Theorem [TT|) that it is possible to hnd a 
solution of the invariance equation in the sense of formal power series (this is a very standard 
order by order perturbation expansion, but in Section [10] we also construct fast algorithms 
that double the order in each step). By truncating appropriately this formal power series, 
we obtain functions that solve the invariance equation very approximately. These approx¬ 
imate solutions are taken as the initial points of an iterative procedure which is shown to 
converge by an a-posteriori theorem (see Theorem IT4ll . which is very similar to Theorem 
20 in jCCdlTdO^ . The domain Q is obtained by examining carefully the process, and the 
quantitative and explicit conditions of Theorem [Ml Notice that for the applications in this 
paper, it is essential that Theorem [TT] is formulated in an a-posteriori format (we can start 
the iteration of an approximate solution even if the problem is not close to integrable and 
also we can obtain estimates on the distance from the initial data to the solution by the 
error of the initial approximation). 

Methods similar to those described above (using a formal power series as a jumping off 
point of an a-posteriori method) were used in |JdlLZ99] . One can also mention |CCdlL13b 
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ICCCdlLTs] which deal with dissipative systems, although in the two latter papers the itera¬ 
tive method is not a Newton’s method, but rather a contraction argument (taking advantage 
of the fact that the dissipation is very strong). The method is very well suited for the study 
of monogenic properties and Whitney extensions, which we study in Section [101 

We expect that the phenomena uncovered here also hold in several other problems. 

1.1. Conformally symplectic systems. As mentioned above, there are many physical 
problems that lead to the study of conformally symplectic systems and, in particular, to the 
study of singular series. 

One simple but very useful remark for conformally symplectic systems is that quasi- 
periodic solutions with enough frequencies satisfy the so-called automatic reducibility: in a 
neighborhood of an (approximately) invariant torus, one can hnd a system of coordinates in 
which the linearization becomes (approximately) constant coefficients. Automatic reducibil¬ 
ity happens irrespective of whether the system is close to integrable or not (see |CCdlLl,^ ). 
Automatic reducibility allows one to develop a KAM theory f |CCdlL13c] ). leading also to 
efficient and accurate algorithms. 

In jCCdlTd3^ one can hnd a numerically accessible method to compute the breakdown 
threshold extending the arguments of [CdlLlOb] . This criterion for breakdown, based on 
the study of the growth of Sobolev norms, is shown to converge to the right value of the 
threshold; indeed, the accuracy of the computed breakdown in actual computers is only 
limited by the available memory, precision and computational time. 

A numerical implementation in actual computers was done in |CC10] . which includes com¬ 
parisons with other methods. Another property of KAM tori in conformally symplectic sys¬ 
tems is that the breakdown of the invariant circles happens due to a bundle collapse scenario 
in which hyperbolicity is lost because the stable bundle becomes close to the tangent bundle, 
even if the exponents remain bounded away (see |CF12j for a numerical implementation and 
a presentation of empirical results, including scaling relations for the breakdown). 

We also mention that the Greene’s criterion for the breakdown of invariant circles has 
been extended to conformally symplectic systems and given a partial justihcation (see 
[GGFdlTddj b 

It is known that KAM theory for conformally symplectic systems requires adjusting param¬ 
eters 1 jCGdlTTScj l and moreover (see [CCdlLTSa]) that Birkhoff invariants near a Lagrangian 
invariant torus with a dynamics conjugated to a rotation disappear (i.e., given a Lagrangian 
invariant torus with a dynamics conjugated to a rotation, there is an analytic and symplectic 
change of variables dehned in a neighborhood of the torus, that conjugates the dynamics to 
a rotation in the angles and a multiplication by a constant in the actions). 

The papers |SL12l ISL15j develop a Kolmogorov theory based on transformation theory for 
quasi-periodic solutions in quasi-integrable conformally symplectic systems and implement 
it numerically. The paper |CC09j develops an a-posteriori KAM theory for the spin-orbit 
problem (a two degrees of freedom conformally symplectic system). 
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1.2. Description of the set up. We will consider analytic families of mappings or flows 
with a small parameter e and having also internal parameters /i. That is, given an analytic 
symplectic manifold M of dimension 2d with exact symplectic form fl, we will consider 
families of mappings '■ M M satisfying: 

= A(£)fi , A(0) = 1 (1.1) 

or families of flows such that: 

, x(0) = 0 , (1.2) 

where A = A(e), y = x(^) fhe conformal factors for maps and flows, respectively, the 
star in fll.ip denotes the pull-back, in fll.2p is the Lie derivative. We will refer to A(e:), 
x{s) as the dissipation. 

Since we will discuss analyticity, all the parameters will be taken to be complex. In 
applications the parameters are often real and then the values of the functions are real. It 
will happen that all the calculations we perform respect the properties that real arguments 
of the function lead to real results. 

In formula (II.ip . e is a small parameter that controls the dissipation. The parameters 
/i G are some intrinsic parameters of the model that are called the drift in some papers. 
Of course, the case when A = 1 (respectively, x = 0) corresponds to the mapping 
(respectively, the flow J>,e) being symplectic. 

Note that the interpretation of (11.11) or (11.21) is that the mapping or the flow transforms 
the symplectic form into a multiple (could be a complex multiple for complex e) of itself. 
We assume that the conformal factor is of the form 

A(e) = 1 + ae^ + 0{\e\^+^) , 

where a G a G C\{0} (in many applications a G M\{0}, so as to preserve the property 
that real arguments lead to real variables). 

We will hx bj within the set of Diophantine vectors (see Dehnition [3] for the standard 
dehnition for maps and Appendix |A] for the standard dehnition for hows). We will be 
interested in studying the domain of complex values of e for which we can continue a KAM 
torus invariant for the symplectic system. 

Of course, we will assume several non-degeneracy conditions (similar to the twist condition) 
to prove the existence of KAM tori. These non-degeneracy conditions can be verihed by 
performing some calculations on the approximate solution. 

Note, that we will not assume that the Hamiltonian system is integrable or nearly- 
integrable, but only that it has a KAM torus of frequency uj. In particular, the results 
apply to perturbations of islands generated by resonances at higher values of the perturba¬ 
tion { |Dua94] . |Dua99] L 

For simplicity we will describe in detail the case of maps, while hows are discussed in 
Appendix 
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The way we seek invariant tori of mappings is to try to find an embedding ^ M. 

and a parameter vector /i^ G in such a way that 



(1.3) 


where is the shift map dehned by a;, T^(6') = 6 + uj. We will hx u to be Diophan- 
tine (see Section ITT]) . Of course, the equation fll.3p will have to be supplemented by some 
normalization conditions, which ensure that the solutions are locally unique. We refer to 
[CdlLOQl ICdlLlOal IdlLRQH ICF12j for a method to hnd invariant curves that solve the in¬ 
variance equation fll.3p numerically. 

Our main result. Theorem [Ill shows that, if there is a solution of (ll.3p for e = 0 (the 
symplectic case), which satishes some mild non-degeneracy conditions, we can hnd and 
/ig dehned for a set Q oi e. The functions and the vectors fis are analytic in e when 
e ranges in the interior of the sets which are described in Section 12.2.11 They also extend 
continuously to the boundary of Q. We anticipate that the sets Q do not include any ball 
centered at the origin in the complex plane, even if they contain the origin in their closure. 
On the other hand, the sets Q fail to include a ball centered at the origin by very little. 
As we will see, the domains are obtained by taking from the ball centered at the origin a 
sequence of smaller balls centered along smooth lines going through the origin and with radii 
decreasing faster than any power of the distance of the center to the origin (see Section 12.2.11 
and Figure HI). 

Similar analyticity domains appeared in [JdlLZ99j . where the authors used a strategy 
close to ours in considering the domains of analyticity of resonant tori in near-integrable 



(see also |GG051 IGGG061IGGGG07] for results based on resummation of series for 


the same problem as |,TdlLZ99] ). The paper |,TdlLZ99j also obtained other geometric results 


such as the monodromy of the stable and unstable bundles, which are not present in other 
treatments. 


We note that the domains we obtain here have several cuts and that the analyticity 
domains do not contain sectors centered at the origin with aperture bigger than vr/a. Hence, 


one cannot use only general complex analysis methods ( [Har49[ IBal941 ISokSOj ) to deal with 


the asymptotic series, as these methods do not guarantee that there is only a function with 
the same asymptotic expansion ( [PLOSl [SZ65] ). On the other hand, it is a byproduct of our 
analysis that the expansions of solutions of the equations using the Lindstedt procedure are 
indeed unique. In Section fl 0.51 we also discuss 0°°-Whitney properties of the solution. 

The argument we present to prove our results has two ingredients: 

{i) An a-posteriori KAM theorem (Theorem [T4P for conformally symplectic systems 
with complex parameters, which shows that if there is an approximate solution of 
the invariance equation fll.3l) . then there is a true solution nearby. 

^The problem considered in |JdlLZ99] is actually more singular than the one considered here, since the 
twist and other non-degeneracy constants depend on e and they degenerate as e —0. Nevertheless, they 
degenerate like a fixed power of e and one can construct approximations to any order. 
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Theorem [H] is a very quantitative statement on when a solution is approximate 
enough to be the starting point of an iterative algorithm which converges. The main 
condition of Theorem [T4l is that the initial error is small enough. The precise smallness 
condition depends mainly on the number theoretic properties of the complex number 
A - the conformal factor - with respect to the frequency uj. The smallness condi¬ 
tion depends also on the Diophantine properties of u and on some non-degeneracy 
conditions of the map. The Theorem [H] is a small modihcation of Theorem 20 in 
[CCdlTd.l^ . 

{ii) An algorithm to produce a perturbative series expansion that provides approximate 
solutions to all orders in e, see part A of Theorem [121 

This algorithm is quite a standard procedure, which goes back at least to [P^ 
and is based on earlier literature. 

In our case, taking advantage of the fact that the maps are conformally symplectic 
(and, hence, automatically reducible) we can improve the classical results in several 
directions: In Section [10.2l we present a quadratic algorithm to compute the Lindstedt 
series, which is much faster and which could replace part A of Theorem [121 Indeed, 
the quadratic procedure in Section [10.21 gives an alternative proof of all of Theorem [T2l 
We also note that using the automatic reducibility we can develop these series starting 
at any invariant torus. This leads to improvements on the domain of analyticity that 
are discussed in Section 110.41 as well as to Whitney regularity in the boundary. 


Theorem [121 is proved by showing that an iterative procedure (which is explicitly de¬ 
scribed in Algorithm [161 and which is a very practical numerical algorithm implemented in 
[CClOl ICF12j l converges if the initial error is small. We point out that, if we start the 
iterative procedure not on one parameterization and a drift, but on an analytic family of 
parameterizations and drifts indexed by a variable e, then the result will also be an analytic 
family in e (the iterative procedure - see Algorithm [161 ^ consists in applying algebraic op¬ 
erations, taking derivatives and solving cohomology equations with constant coefficients; all 
these elementary steps preserve the analytic dependence on parameters). We also point out 
that the convergence will be uniform in a domain of £, if all the non-degeneracy conditions 
and the initial error are uniform. 

Hence, if we start the iterative procedure with a function analytic in an open set (respec¬ 
tively, continuous in a closed set) of e (such as that produced by part A in Theorem [T^ . 
then we obtain that the limit is analytic (respectively, continuous) in e in the domain where 
the convergence is uniform. 

Since, furthermore, we have local uniqueness of the normalized solutions of the invariance 
equation that satisfy a normalization condition, we are sure that the solutions obtained in 
two open sets agree on the overlap, and hence we can use arguments based on analytic 
continuation. 
















R. CALLEJA, A. CELLETTI, AND R. DE LA LEAVE 


For simplicity, we will present the proof for maps and we provide the changes needed to 
obtain the result for flows in Appendix 

This paper is organized as follows. In Section |2l we collect some of the standard dehnitions; 
in Section 12.11 we provide some dehnitions on Diophantine properties, while in Section 12.21 
we present the geometry of several sets where the conformal factors satisfy Diophantine 
conditions. In Section [3] we present the invariance equation and the normalization conditions 
we impose to obtain local uniqueness. 

In Section 0] we state Theorem [121 fhe main result of this paper. 

In Section[5]we present Theorem [T4l which is the a-posteriori theorem alluded in (i) above. 
Such a theorem is the hrst ingredient of the main result. The proof of Theorem [TT] is very 
similar to the proof of Theorem 20 in |CCdlL13^ and we just outline the (rather minor) 
differences. The existence of the Lindstedt series is discussed in Section [6l In contrast to 
the procedure using Lindstedt series f |CCdlLl^ L the present treatment allows one that the 
conformal factor depends on e. This provides the proof of the hrst part of Theorem [121 while 
the second part is proved in Section [3 In Section [8] we study some geometric properties of 
the set. In Section [9] we include a conjecture on the optimality of the results described in this 
work. In Section [10] we present several results: we show that the Lindstedt series expansions 
can be obtained around any point; we present a relation with the theory of monogenic 
functions and establish that the embedding function and the drift are Whitney differentiable 
in the domain; we provide a quadratic method for the computation of the Lindstedt series, 
leading also to Part B of Theorem [121 we present a discussion on the improvement of the 
domain of analyticity. The extension to the case of flows is provided in Appendix [Al 

Remark 1. In many estimates in this paper, we can obtain that the domains satisfy upper 
bounds less or equal than for all N and for some constant C^. Clearly, since the 

bounds are true for all N, one can get an upper bound less or equal than inf at = r(e). 

The function T will, of course, go to zero faster than any power. 

If one had explicit forms for Cn, it would be possible to obtain explicit forms for T. In 
many problems similar to the ones we are considering, one obtains factorial bounds like = 
AN°‘^ for some constants a and A. In such a case, one obtains that r(£) = exp(— 
for some constant C. We conjecture that the sizes of the balls to be excluded in the present 
paper satisfy the factorial bounds. 


2. Some definitions 

In this Section, we collect some dehnitions on spaces, Diophantine properties and we 
set the notation. Most of the dehnitions in this section are standard. One non-standard 
dehnition that will play an important role is the Diophantine property of complex numbers 
with respect to a Diophantine frequency (see Dehnition 01). 

Given p > 0 we dehne the complex extension of the d-dimensional torus as 

T^p = {ze : Re{zj) G T , | lm{zj)\ <p , j = 1, ..., d} . 
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We define Ap to be the vector space of functions analytic in Int(T^) and which extend 
continuously to the boundary of T^. 

We endow Ap with the supremum norm 


= sup \f{e)\ 


e&Tf 


( 2 . 1 ) 


The norm fl2.ip makes the space Ap into a Banach space, indeed a Banach algebra under 
multiplication. An important closed subspace of Ap is the set of functions which take real 
values for real arguments. 

Analogous dehnitions are made for analytic functions on taking values in vectors or in 
matrices and, of course, the multiplicative properties of norms of vectors and matrices lift 
to supremum norms of functions taking values in vectors or matrices. 

We also note that we can dehne analytic functions of e G C taking values in Ap. Following 
standard dehnitions, we say that an ^p-valued function K is analytic in an open domain 
P C C when, for any Eq G V, we can represent for all |e — £o| sufficiently small, = 
, where the convergence of the inhnite sum happens in Ap. Note that, with 
this dehnition, it is clear that an ^p-valued analytic function is also an ^p/-valued analytic 
function for p' < p. 

It is remarkable that there are many other dehnitions that are apparently weaker than 
the dehnition above, but which turn out to be equivalent (see |Hil48l Chapter III]). In our 
case, the K will be analytic functions from Q to some Ap. In principle, the domain Q could 
depend on p, but we will not include it in the notation. 

We note that we are only obtaining lower bounds of the domain of analyticity. The method 
of proof leads to estimates for diherent p’s, which are proved by selecting a diherent constant. 
Hence, the statements that are valid for all Diophantine constants in a range are also valid 
for all the values of p in a range. 

We recall the classical Cauchy inequalities for derivatives and for Fourier coefficients (see, 
e.g., |SZ651 lRhs75] ). 


Proposition 2. For any 0 < 5 < p and for any function K G Ap, denoting by the j-th 
derivative, we have: 

\\D^K\\p_s < C,5-^\\K\\p, 

\Ku\ < ||iF||pe-2-l"l^ 

for some constants Cj, where denotes the k-th Fourier coefficient of K and \k\ = |fci| + 
... + \k(i\. 


2.1. Diophantine properties. In this Section we collect some dehnitions concerning Dio¬ 
phantine properties, which will be needed to bound the small divisors appearing in the 
solution of the invariance equation. The only non standard dehnition is Dehnition 01 

Definition 3. Let u G r G M+. We define the quantity u{u}]t) as 

u{u;t)= sup — l\~^\k\~^ . (2.2) 

fcez<*\{o} 
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In the sup above we allow oo. Also, if — 1 \ = Q, we set z/(a;; r) = oo. 

We say that cj is Diophantine of class r and constant v{uj]t), whenever 

iy{uj-, t) < oo . (2.3) 

We denote by Prf(z/, r) the set of Diophantine vectors in of class r and constant p. 

Of course, if u is Diophantine of class r, it will be Diophantine of class ^|J for all > t. For 
the purposes of this paper, the value of the constant u is more important than the exponent 
r, so we will consider one hxed exponent in the main theorems. 

Definition 4. Let A E C, cj E r G M+. We define the quantity z/(A; cn, r) as 

p[\-,u:,t)= sup . (2.4) 

fceZ'^Uo} 

Again, we allow oo in the supremum above and set p{X;u,t) = oo, if = A for some 

k E 

Remark 5. We note that, for a fixed oj, the function v = is a lower semi- 

continuous function of X, since it is the supremum of continuous functions. 

Remark 6. If z/(A; u, r) < oo, then for all > r, we have i/(A; < oo. 

Note that the definition of p{X;u), r) does not require that to is Diophantine of exponent r. 

Remark 7. No matter what u is, if \X\ ^ 1 , then from the inequality 

_ ^||^|r > |g2«fc.a; _ > 

one obtains that 

p{X]u,t) < |1 — |A||“^ < oo 

for any r. 

As a consequence of the above remark, the only case that needs to be studied in detail 
is when |A| = 1, in which case it could be that p{X-,u,t) = oo (e.g., if A = for some 

k E'L'^\ {0} or if cj is a Liouville number). 

For T > d — 1 the set of real vectors which are not Diophantine of class r has zero Lebesgue 
measure in M'^. That is, the union over all z/ > 0 of the sets of Diophantine vectors of class 
r satisfying (12.31) with constant z/, has full Lebesgue measure in M'^. For any uo E M'^, if 
r > d — 1, the set of A for which z/(A;a;, r) < oo is of full Lebesgue measure on the unit circle 
(see |Sch8nj ). When we consider complex to E VA, the set of vectors that are not Diophantine 
of class r has zero measure when 2r > d — 1. This shows that complex Diophantine is easier 
than Diophantine over the reals. 

We will consider hxed the analytic function A(e), which gives the conformal factor as a 
function of the perturbing parameter £. In particular, we will assume that A(£) is analytic 
in a neighborhood of zero and that A(0) = 1. Hence, we assume that A = A(e) satishes: 


^2'Kik-U) I 


-|A| = 1-|A| 


HA 


A(e) — 1 = ae°‘ + 0(|£|“+^) 
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for some a > 0 integer, a G C\{0}. 

This will be one of the assumptions of our main result stated in Theorem IT^ Note that, 
since we are considering analytic functions, the alternative to the existence of a and a in 
HA is that X{e) = 1, that is that the maps are symplectic. Hence, given the analyticity 
assumption, the only content of HA is that the perturbations indeed change the symplectic 
character (as well as setting the notation for a, a, which will play a role in the the quantitative 
statements). 

We consider the function A(e) as given in HA. 

We dehne the set Q = GiA; u, r, N) for some H > 0, iV G a; G as 

G{A;u),t, N) = {e e C : i^(A(£);a;, r)|A(e) — < H} . (2.5) 

We also introduce the notation 

Gro{A-,uj,r,N) = G n{e e C |e| < ro} ; (2.6) 

the set Gro will typically be used for sufficiently small vq. 

Notice that G{A]U}, r, N) is the preimage under the function A of the set 

A(H; w, r, AT) = {A G C : z/(A;a;, r)|A - < H} . (2.7) 

The motivation for introducing these sets will be discussed in Section 12.21 We anticipate 
that the set G is the set where the Diophantine constants of A(e) are not too bad, so that a 
good approximation (up to a high power of e) can be taken as the initial condition for an it¬ 
erative procedure that converges. More detailed motivations will be presented in Section 12.2( 
some depiction of the sets Gro: A is presented in Figured! 

Since, for a fixed u, the function u is lower semi-continuous, we note that ^ is a closed 
set. The interior of this set is non-empty and indeed, the set G is the closure of its interior. 
We will be studying functions dehned in G taking values either in the complex or in some 
Banach space of functions. We will consider functions on G which are continuous in G and 
analytic in the interior of G- We will sometimes refer to these functions as analytic in G- 

For obvious typographical reasons, since many of the arguments of G, A will be hxed in 
the discussion, we will omit them: if u, r, N, A are hxed in an argument, we will just write 
GiA) or Gro{A). Of course, one could also add p to the dependence of G, but as remarked 
above, this is not needed. 

2.2. Motivation for the role of the Diophantine constants. In this paper, we will 
hx uj Diophantine and search for solutions of the functional equation fll.Sp . capturing the 
geometric idea that we have a parameterization of an invariant torus. 

The solutions of the invariance equation fll.Sp are obtained by an iterative method, whose 
step involves solving two cohomology equations of the type considered in Section l2. 2.2l below 
(as well as algebraic and calculus operations). 

One of the cohomology equations in the iterative step will involve small divisors of the 
form ~ appearing in (12. 2 p and the other will involve small divisors of the form 

^^ 2 mk-uj _ ^|-i appearing in fl2.4p . Hence, the quantitative hgure of merit of the step will 
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be the constant r)z/(A; cj, r) for some r, f G M+ (as well as other factors, that we will 
consider hxed). 

Since the two cohomological equations that we are going to solve in Section |5]are different, 
there is no reason to impose that the exponents are equal, but, for the sake of simplicity, we 
will not optimize these choices and, from now on, we will just consider the case f = r. 

Based on the results presented in |CCdlLl,^ . we will show that, if we start with an 
approximate solution of fll.Sp and if the initial error is small enough compared with the 
hgure of merit (the precise relation given in assumption H5 of Theorem IT^ . then the iterative 
procedure can be repeated inhnitely many times and it converges to a true solution. Taking 
as initial condition the Lindstedt series, we will obtain convergence in the sets claimed in 
Theorem [I2l 

If we £x the frequency ca and the exponent r, the quality factor for the step, namely 
1 /( 0 ;; r)z/(A; (V, r), is a function of A alone. It will be important to study the places where 
this function is not too big and identify complex domains in the £-plane where it is bounded 
uniformly. The values of A which lead to a quality factor in the iterative step, which is small 
enough with respect to the initial error, are those for which the procedure converges. 

As indicated in hypothesis HA, when we consider problems depending on the complex 
parameter £, the quantity A(e:) — 1 will be approximately «£“. The perturbation expansion 
to order N will produce approximate solutions that satisfy the equation to order (see 
(14. 2 p in Theorem [12]) . Hence, the sets Q in (12.5p are the sets where we can start the iterative 
procedure with an approximate solution given by a perturbation expansion, and obtain 
that the iteration process converges. The sets Q in (12.5p will thus be the sets where we 
can establish that there is a solution starting the iterative procedure from the approximate 
solution given by the expansion. Of course, it is possible that the set of e for which K^, pg 
are analytic is larger than Q. The set Q is the place where one particular method works, but 
we could use other methods. Indeed in Section 110.41 we will find some way to extend the 
domain Q. 

We will undertake the study of the geometry of the sets Q and A in Section 12.2.11 and in 
Section jH] We anticipate that, of course, one can obtain uniform bounds for z/(A; uj, r) in 
any domain which is away from the unit circle. The function u{X]u, r) is inhnite in a dense 
set on the unit circle, but can be bounded in domains that include the unit circle in the 
boundary, indeed in a full measure set in the unit disk. The geometry of these domains is 
rather interesting and will be discussed in the next Section. 


2.2.1. Geometry of the level sets of u{X;u),t) . In this Section we try to understand the sets 
of points described in Section 12.11 More properties will be described in Section |H1 

The complement of the set A, say A'^, is the bad set where the bounds we use do not allow 
us to claim the existence of tori, but we claim that if we take as initial condition of the 
iterative procedure the Lindstedt polynomial of order N, then the measure of A'^ is small for 
N large enough. 
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The set A'^ is the union of the open sets 

7^fc = {A e C : > A\kY\\- l|-(^+h} 

= {A e C : - A| < A-^\k\-^\\- 1|^+^} . 

Since A appears on both sides of the inequality dehning IZk, this is not easy to handle. On 
the other hand, we note that if we consider the annulus p < |A — 1| < 2p, we see that the 
intersection of the set TZk with the annulus is contained in the ball Bk with center at 
and radius CA~^p^^^\k\~'^, for some constant O > 0. 

For each k {0}, the area Sk of the ball is equal to 

Sk = 71 A-^ p2(A+l) |fc|-2r _ 

Notice that the area of the union of the balls Bk is hnite, provided 

\k\~^^<oo. 

k£Z‘^\{0} 

Indeed, we see that as we consider p small, the excluded area decreases faster than 
which is much smaller than the area of the annulus, which is proportional to p^. Hence, we 
can say that A has a point of density at 1 and that the excluded balls in the A-space decrease 
very fast as we approach 1. More detailed estimates will appear later. 

The above remarks can be translated to the e plane under the assumption HA. 

Figure [U provides a representation of the circles which must be excluded from the dehnition 
of the sets introduced in (12. 5 h and (12.7p . We remark that taking a higher exponent N in 
fl2.5p . fl2.7p does not alter signihcantly the geometry of the domains, it just makes the radius 
of the circles decrease faster. This fast convergence to zero of the radii makes it impossible 
to represent the excluded regions in a quantitatively correct way. All the circles except for 
a few of them will be smaller than a pixel! Nevertheless, we conjecture that these excluded 
circles are there, see Section m 

2.2.2. Estimates of cohomology equations. In this Section, we will present (rather elemen¬ 
tary) estimates on the solutions (p : T'^ —)■ C of twisted cohomology equations of the form: 

Xip{9) - ip{9 + u) = p{9) , (2.8) 

where the function r; ; T'^ —)■ C, the parameter A G C and the frequency vector ca G are 
given. 

We will assume that oo is Diophantine (see Dehnition [3]) and that A satishes Dehnition |4] 
with i/(A; uj, r) < oo. We want to show that there exist solutions of fl2.8p and that we obtain 
quantitative estimates on the size of these solutions in terms of the quantitative estimates of 
the Diophantine properties of A. The estimates that we will obtain will be tame estimates 
in the sense of Nash-Moser implicit function theorems. 

For subsequent applications in this paper, it will be quite important that the estimates 
that we obtain about (f are rather explicit on Hence, they will hold uniformly in 

sets of A for which z/(A; u, r) is uniformly bounded. The geometry of these sets was described 
explicitly in Section 12.2.11 
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Figure 1. A representation of the sets A, Q introduced in (12.Th . (I2.5p : they 
are given by the region not covered by the black circles. The domains are 
obtained for d = 1, t = 1, a = 3. The radii of the balls have been rescaled 
for graphical reasons. Left; in black the complement in the complex plane of 
the set A as in fl2.7p ; the centers of the balls he approximately on the unit 
circle, which at high magnihcation resembles a vertical line. Right: in black 
the complement of the set Q as in fl2.5p : the centers of the balls are at “ — 

i|i 

Since we will be interested in the limit A ~ 1, it will be important that we can consider A 
ranging on domains which include 1 as a limit point. Hence, we will insert in our hypotheses 
that uj is Diophantine to avoid empty statements. The following Lemma is a standard result 
in KAM theory; the dependence of the estimates on the exponents is not optimized as it is 
done for A real according to |Riis751 [RiisTGaj . This will lead to slightly different estimates 
in H5 of Theorem [T3] when compared to those appearing in Theorem 20 of |CCdlLl,^ ; a 
more detailed comparison will be given later in Section [5l 

Lemma 8. Let u G t), A G C. Assume that 7] G Ap, p > 0, is such that Ti{0)d9 = 0. 
Then, we can find a unique p G L^(T^) solving fl2.8l) . that also satisfies 



Moreover, for any 0 < 6 < p, we have p G Aps, and 

\\<p\\t-s < C{T,d) V"'* IIijIIp 


(2.9) 
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for a suitable constant C = C{T,d). 

Note that a particular case of Lemma |8] is when A = 1, which reduces to the most standard 
cohomology equations of KAM theory. 

Proof. The proof we present here is based on the most elementary (but possibly not optimal 
in the exponent for 6) argument (see Remark [9] below). 

We expand 7] in Fourier series as ri{9) = obtain that fl2.8p is 

equivalent to having that for all k & if \ {0}: 

X,-! ^27Tik‘UJ,'x _ ^ 

^ ^k Vk 5 

whose solution is 

fk = {X- . 

We then have the following estimate: 


p—S — 




k ^ 


27T{p—5)\k\ 


< 


E 1 ^ 


^2'Kik-ijj 1—1 


\flk\e 


2'K{p—5)\k\ 


fcGZ^\{0} 

< ^ z/(A;u;,r)|/c|’ 

feGZ‘*\{ 0 } 

fceZ'^\{o} 

< Cn{\-uj,T)\\7]\\p 


k&d\{0} 


^-2-Kp\k\ 2-K(p-6)\k\ 


3 —27r5|fc| 


E^' 


r+d—l^—2'K6j 


where we have just used the Cauchy estimates for the Fourier coefficients in terms of the 
supremum in a band and the dehnition of the constant z/(A;a;, r) in fl2.4p . The desired result 
(12. 9 p just follows from estimating the last sum, which is easily shown to be asymptotically 
bounded by □ 


Remark 9. When A G M, the papers |Rus75[ IRiisTGb] contain more sophisticated esti¬ 
mates that lead to conclusions in Lemma 0 with a better exponent on 6. Namely, with 
the same notation of Lemma 0 when A G M, |Rus751 IRus76b] lead to the conclusion that 
|| 7 ^||p -5 E Ck'5~'^\\r]\\p. Note that not only the exponent in |R,iis75] is better than the exponent 
in Lemma 0 but also that the constant appearing in the bounds above is proportional to the 
Diophantine constant of u. 

Unfortunately, when A G C \ M, it seems that it would be necessary to reexamine carefully 
the proof of |Rus75l IRus76a] . This would indeed be quite interesting on its own merit, but 
will not be pursued here since, for the purposes of this paper, the geometry of the level sets 
of the constants in the estimates in Lemma 0 is more important than the exponents. The 
straightforward argument above leads to constants n, whose geometry is easy to analyze. 

Note that the Diophantine properties ofu enter into the proof of Theorem 20 in |CCdlL13c] 
only through the estimates of the linearized equations. Hence, using the elementary estimates 
in Lemma 0 is equivalent - for the proof of Theorem 20 in |CCdlLl,^ - to considering u 
with a different Diophantine exponent. 
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Remark 10. We recall that, due to Remark^ when |A| 7 ^ 1, there are no small divi¬ 
sors in the solution of the cohomology equation fl2.8p . Hence, when |A| 7 ^ 1, we have that 
Mlp-s < C'(||A| - l|)-i 6-^-^ \\r]\\p. Indeed, ||(^|lp -5 < C{\\X\ - 1|)-^ \\r]\\p. 


3. The invariance equation and some normalizations 


3.1. The invariance equation. The centerpiece of our treatment will be the invariance 
equation 

fp^,,oK, = IQoT^ . (3.1) 

We think of (13.ip as an equation for the parameters /ig and for the embedding once £ is 
hxed. Of course, when e varies, we can think of K, ji as functions of e. 

We will develop a Newton’s method that can start from any approximate solution. Using 
the geometric properties of the system, the equations can be reduced to constant coefficients. 
This has been the idea of the KAM theory based on parameterization, started in |dlLnip 
IdlLGJVO^ for symplectic systems; the extension of the formalism for conformally symplectic 
systems is developed in |CCdlL13^ . 


3.2. Some normalizations and uniqueness. We note that the equation fl3.ip never has 
unique solutions. If {K,ia) is a solution so is, for any a E T'^, (A'^°'\/i) where = 

K{9 + cr). Of course, describes the same torus as K, the only difference is the origin of 
the parameterization. So, we can hope to get uniqueness only if we impose a normalization 
that hxes the origin in the 6 variables. Later we will see that, indeed, this change of the 
origin is the only source of non-uniqueness and that, once we £x it in any reasonable way, 
we obtain that the solutions of (13.11) and the normalization are locally unique. For us the 
normalization is important because it allows for analytic continuation of solutions produced 
by different methods. 

The following normalization has been found to be geometrically natural and easy to im¬ 
plement numerically. 

Given a reference embedding Kq, which we choose once and for all, we can form the matrix 
M obtained juxtaposing the matrices DKq and J~^ o KqDKqN, say 

M = [DKo I J-^ o KoDKoN] , (3.2) 


where 

N = {DKj DKo)-^ (3.3) 

is a normalizing matrix. As we will see below (see |GGdlLl,^ for more details), the matrix M 
provides a frame of reference near the image of Kq. Notice that DKq transforms vectors that 
are in the tangent space to the image of Kq and that, since Kq{T‘^) will be almost Lagrangian, 
we have that J~^ o KqDKq - the symplectic conjugate - will be almost perpendicular. 

We will say that a torus with embedding Kq is normalized, when 


Ifd L 


M-\e){K{e) - Kq^o)) 


de = t), 


(3.4) 
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where the subscript 1 denotes that we take the hrst d rows. 

The geometric interpretation is that M dehnes a particularly interesting system of coor¬ 
dinates near the torus Kq. The quantity M~^{6){K{6) — Kq{6)) expresses the displacement 
from Kq{6) to K{0) in this system of coordinates and our condition is that, in this system 
of coordinates, the 6 component of the displacement has zero average. 

With reference to the normalization (I3.4p . we recall also the following result (Proposition 
26 from |CCdlL13^ ). which shows that it can be imposed without any loss of generality for 
solutions that are close. 

Proposition 11. Let Ki,K 2 be solutions of fl3.ip . \\Ki — iP 2 |lci be sufficiently small (with 
respect to quantities depending only on M - computed out of Ki - and f). Then, there exists 
a e M"-, such that = K 2 oT„ satisfies fl3.4p . Furthermore: 

|c’‘| < C* ||Ai — iP2||co ) 

where the constant C can be chosen to be as close to 1 as desired by assuming that Ki, 
K 2 are twice differentiable, DKfDKi is invertible and \\Ki — /P 2 IIC 0 is sufficiently small. 
The a thus chosen is locally unique. 

Note that Proposition fTTl shows that if we have a family of solutions close to Kq, we can 
modify this family of solutions by composing them with a small displacement, so that they 
are normalized solutions. 

4. Statement of the main result. Theorem [12] 

In this Section, we state the main result. Theorem [121 which - under suitable assumptions 
on the mapping, the freqnency and some non-degeneracy conditions - allows us to prove the 
existence of an exact solution of the invariance equation, analytic in the set introdnced 
in ([23]). 

Theorem 12. Let M. xB with B an open, simply connected domain with smooth 
boundary; M is endowed with an analytic symplectic form 12. Let us denote by J = J{x) the 
matrix representing 12 at x, so that for any vectors u, v, one has 

flx{u,v) = {u,J{x)v) . 

Let cj G satisfy Definition 0 for Diophantine vectors. 

Let f^^s with /i G P, P C open, e E C, be a family of conformally symplectic mappings, 
that satisfy fll.ll) with conformal factor as in HA. 

Assume that for e = 0 the family of maps ffj,^ is symplectic. 

Assume that for some value Ho the map f^^^ admits a Lagrangian invariant torus, namely 
we can find an analytic embedding from —)■ Ai, say Kq G Ap(T‘^, Ai), such that 

f^J,o,o ° Kq = Kq o 3 ^ , 

where Kq G Ap for some p > 0. Moreover, assume that Kq is Lagrangian, namely that 

DKj J o Kq DKq = 0 . 


(4.1) 
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Assume furthermore that the torus Kq satisfies the following hypothesis. 
HND Let the following non-degeneracy condition he satisfied: 


det 


So Sq{Bi ,q)^ + Aiq 
0 A20 


7^ 0 , 


where the d x d matrix So is defined as 


So{e) = N{e + ujy DKoie + uj) o K^ie) j-^ o K,{e) DKo{e)N{e) 

- N{e + uYDKo{e + uYj-^ o Ko{e + u) dko{o + uj)N{e + u) 

with N as in (EJl), the d x d matrices Aw, I 20 denot^ the first d and the last d rows of 
the 2d X d matrix Aq = (M o (-D^/mo.o o Kq), where M is as in fl3.2p . (i?bo)° is the 

solution (with zero average) of the cohomology eguation (Bbo)^ — (Bbo)^ = —(^ 20 )°, where 
{BboY = Bho ~ Bbo and the overline denotes the average. 

Then, we have the following results. 

A) We can find a formal power series expansion satisfying fl3.ip in the 

sense of formal power series. 

More precisely, defining = Yl!j=o^^h-i Sf eN and p > 0, 

we have 

II f [<.v, ^ O O T^\l, < CN\e\^+^ (4.2) 

f-^E 5^ 

for some 0 < p' < p and > 0. 


B) We can find a set Qro of the form introduced in fl2.6l) with vq sufficiently small and for 
any 0 < p' < p, we can find the functions 

K : Qro -^p' ) 

P '■ Gro , 

which are analytic in the interior of G^o and extend continuously to the boundary ofGro, 
that for e E Gro they satisfy exactly the invariance eguation 

Ue,eoK,-K,oT^ = 0 . (4.3) 


Furthermore, we have that the solutions thus found have the formal series provided in part 
A) as an asymptotic expansion. That is, for any iV G N and for any 0 < p' < p: 


||/cpi-iC,||p < CM\e\^+\ 



(4.4) 


^We call attention that in |CCdlL13c] the statement of Theorem 20 defines the matrices Ai, A 2 as the 
first d and the last d columns of the 2d x d matrix A. Clearly, this sentence does not make sense, unless 
one changes “columns” by rows. Since the sentence as written in |CCdlL13c] is clearly impossible and the 
detailed calculations are given, we hope that this has not misled the readers, but we take the opportunity 
to set the record straight. See also the discussion in Section O 
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Remark 13. The condition HND has a very transparent geometric interpretation, which 
we will present in the proof given in Section\^ See also |CCdlL13^ for more details. 

We also remark that the formal power series can be chosen to he normalized with respect 
to Kq. The can also he chosen to he normalized. 

5. Quantitative a-posteriori KAM theorem for conformally symplectic 

SYSTEMS 

The goal of this Section is to state a quantitative KAM theorem, namely Theorem [Ml 
which is very similar to Theorem 20 in |CCdlL13(^ . We will also detail the (rather minimal) 
changes to the proof of Theorem 20 in |CCdlL13^ . needed to obtain a proof of Theorem [Ml 

Theorem 20 in |CCdlL13(^ shows that, given a fixed family of conformally symplectic 
mappings and an approximate invariant torus for a value of the parameter, we can find an 
exact invariant torus for a nearby value of the parameter. The result presented in |CCdlL13^ 
is based on an a-posteriori format, which is very natural, because we do not need to start 
from an integrable system. 

The main novelty here with respect to Theorem 20 in |CCdlL13c] is that we highlight the 
dependence of the results on A - the conformal factor - and that we allow this conformal 
factor to be complex. 

In |CCdlL13(^ the conformal factor A was considered essentially fixed. In some results of 
|CCdlL13^ . A was allowed to range over a real interval [1 — A, 1 -|- A] for some A > 0. In this 
paper, however, A changes and the main goal of this Section is to consider the dependence 
on A. 

The proof of Theorem [Mj will be just walking through the proof in |CCdlLl,^ . but keeping 
track of the dependence of the constants on A. We remark, however, that the dependence in 
A comes only through the Diophantine constant u{X]u,t) as in Definition [H 

The present treatment requires only minor differences with that of |CCdlLl,^ . More 
precisely, 

• In [CCdlLl,^ there is a treatment both of the analytic and the hnitely differentiable 
case. In this paper, we will only present the analytic case, since it is the one used in 
the applications of this paper. 

• In |CCdlLl,^ there is a treatment both of the case uniform in A and the case for 
hxed A. Moreover, in |CCdlL13^ the parameter A was supposed to be real. In our 
case, we will allow A to be complex. We will obtain estimates that, in the language of 
[CCdlLl^ . are uniform in A, but we will pay attention to the Diophantine constants 
and the geometry of the sets in the complex where these constants take values. 

• In this paper we are using only crude bounds for cohomology equations, whereas in 
|CCdlL13^ we used the Riissmann estimates |Rus75] . We do not know how to adapt 
the Riissmann estimates to the case that A is complex (see Remark [9]). 

Theorem 14. Let A1 = x B with B an open, simply connected domain with smooth 
boundary, endowed with a scalar product and a symplectic form D. 


































20 


R. CALLEJA, A. CELLETTI, AND R. DE LA LEAVE 


Assume that the following hypotheses H1-H2-H3-H4-H5 are satisfied. 

HI Let cu G 6e Diophantine of class r G M+ and constant u{uj] r). For A G C assume 
< oo, where is defined in ( 12 . 41 ) . 

H2 Let with // G F, F C C'^ open, e E C, be a family of (complex) conformally 
symplectic mappings with respect to a symplectic form 12, that is = A(£)f2 (see flLip ) 
with A(£) complex. 

Let Ka : —)■ /ia £ be an approximate solution of fll.Sp . such that 

fiia,e°Ka-KaOT^ = E (5.1) 


with error term E. 

H3 Assume that the following non-degeneracy condition holds: 


det 


(A - 1) Id I 2 


where the d x d matrix S is defined as 


7 ^ 0 , 


(5.2) 


^(0) = N{e + ufDKfie + cu) o Kfie) J-I o Kfie) DKfie)N{e) 

- X N{e+ ujyDKa{0+ ujyj~^ o Ka{e+ UJ) DKa{e+ uj)N{e+ UJ) (5.3) 

with N as in fl3.3p with Ka replacing Kq, the d x d matrices Ai, A 2 denote the first d and 
the last d rows of the 2d x d matrix A = (M o Ta,)~^ ° Ka), where M is as in fl3.2l) 

with Ka replacing Kq, is the solution (with zero average) of the cohomology eguation 

X{Bb)^ — {Bb)^ oTa, = — (^ 2 )°, where (Bb)^ = Bb — Bb and the overline denotes the average. 
We denote by T the guantity 

S 

(A - 1) Id A2 

and we refer to T as the twist constant. 

Assume that Ka G Ap for some p > 0. Assume furthermore that for p, eT we have that 
fp^e is a -family of analytic functions on a domain - open connected set - C C LA/ifi x LA 
with the following assumption on the domain. 

H4 There exists ( > 0, so that 




dist(/ia,9F) > C , 

dist(iF,(T''),aC) >C . 


Finally, let the solution be sufficiently approximate according to the following assumption. 
H5 For some 0 < 5 < p, the error term E in fl5.ll) satisfies the ineguality 

2 


\p<C 


z/(a;;r)z/(A;a;,r) , 


where C denotes a constant that can depend on r, d, T, ||Z2iFa||p, |liV||p, \\M\\p \\M ^||p as 
well as on ( entering in H4. 
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Then, there exists /i^, such that 

o K^ - K^oT^ = t) . 

The quantities K^, fj,^ satisfy the inequalities 
\\K,-Ka\\p-8 < 

l/^e l^a\ ^ 11-^| Ip 5 

for positive constants Ck, C^. 

For a fixed value of £, say e = e*, the solutions of fll.Sp are locally unique, provided the 
normalization condition fl3.4p is satisfied. Following |CCdlL13^ . we have the following result. 

Lemma 15. Let be solutions of fll.Sp with sufficiently 

small, close enough. Assume that satisfies fl3.4p . that the non-degeneracy con¬ 
dition H3 is satisfied at Qjid fj^at the assumption H4 on the domain is satisfied. 

Then, there exists a G such that iF^) = iF^) oTg., pd) = pd)^ 

5.1. Some notes on the proof of Theorem 1141 and Lemma 1151 Theorem [T^ is proved 
in |CCdlL13(^ . but without keeping track of the difference between z/(ci;;r) and z/(A;a;,r), 
since only a fixed A was considered in |CCdlL13^ . 

For the sake of completeness, we will just repeat the main steps, so that we can trace 
the constants and verify that the constants are those we claimed in Theorem [TH We omit 
details that can be found in |CCdlL13(^ . 

The procedure of |CCdlL137^ is just to device and estimate a Newton-like method based 
on some identities obtained from the geometric properties of the map. Using these geometric 
identities, the equations appearing in Newton’s method are reduced to constant coefficient 
equations of the form appearing in Lemma [H This procedure is sometimes called automatic 
reducibility. Beside leading to a proof of the reducibility, it leads to an efficient numerical 
algorithm that was implemented in |CC10l ICF12] . An important feature of the method, 
crucial for our purposes, is that it can start from an approximate solution, even if the map 
/ is far from integrable. 

The Newton’s method for equation fl3.ip is based on the following steps. Given an ap¬ 
proximate solution {Ka, La) of the invariance equation (ll.dh as in (15.ip with error term E, 
then find corrections A and a to iF^, /!„, respectively, in such a way that: 

o iF,) a - a o oK,y = -E . (5.4) 

Then, it can be proved that Ka -|- A, /i^ -|- a will satisfy the invariance equation with a much 
better accuracy (in a slightly smaller domain), precisely one expects that the new error will 
be controlled by the square of the old one. 

Unfortunately, the equation fl5.4p is not so easy to deal with because it has non-constant 
coefficients. The idea of the automatic reducibility method is to find an adapted frame of 
coordinates that takes advantage of the geometry of the system. Let us introduce M like in 
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fl3.2p with Kq replaced by Ka'. 

M = [DKa I J-^ o KaDKaN] 

with J the symplectic matrix and N the normalization factor dehned like in fl3.3p : 

N= [DKjDKaY^ 

The geometry of M is that it is a change of basis in the tangent space of the approximately 
invariant torus. The remarkable thing is that, in this basis, o has a very simple 

expression, namely 

Df,^,,oKaM = MoT^(^ 

where R, measuring the error of the automatic reducibility for an approximately invariant 
torus, is a quantity that can be estimated by E (in the sense of Nash-Moser, we allow that 
the estimates are tame estimates in a smaller domain of analyticity). 

The geometric reason for the identity flS.Sp is that, if we take derivatives of flS.ip . we 
obtain the hrst d columns of (15.bp . Geometrically, we have found a vector held DKa that 
gets mapped into itself by the transformation o Ka- The directions J~^ o KaDKaN 

are the symplectic conjugates to DKa (here one uses the fact that the torus is approximately 
Lagrangian, which is established as a consequence that Ka^^a satisfy (I5.ip ). 

Using fl5.5p and writing the correction A for Ka as A = MW for some function W, we 
see that, ignoring a term containing the factor RW, equation (15.dp becomes 

( ^^^W-WoT^ + {MoTa,)-^D^f,^^,oKaa = -{MoTg,)-^E . (5.6) 

Note that fl5.6p is an equation for both W and a. Similar equations appear in KAM theory 
all the time. Note that (15. 6 p becomes a cohomology equation of the form (12. 8 p for the second 
component IU 2 of W and, once this is solved, we substitute W 2 in the equation for the hrst 
component VUi, which is an equation of the form fl2.8p with A = 1. 

Writing (15.6p in components (i.e. taking the hrst d rows and the last d rows) we obtain 

fUi-WioT^ = -SW 2 -E 1 - 
XW2-W2oT^ = -E2-A2a, (5.7) 

where E = (M o Taj)~^E, A = (M o o Ka). We write A as A= [A 1 IA 2 ], where 

Ai, A 2 denote the hrst d and the last d rows of the 2d x d matrix A. 

The solution of equations (15.71) requires that the average of the right hand side is zero, 
but this can be accomplished by properly choosing the quantity a, provided that the non¬ 
degeneracy condition (15.21) is satished. 

This procedure is very standard in KAM theory, but in this case it has a complication. 
If we change a, since it is multiplied by a non-constant function, we change the solutions 
of the equations which we have to seek for the average W 2 and hence (compare with Step 
9 in Algorithm fTbl below) we change the solution for the average in the equation for Wi- 
Therefore, the equations for IU 2 and a are not completely decoupled. The observation in 
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jCCdlLl.Sc] is that we know that the dependence of W 2 on a is affine and that we can 
compute the coefficients by solving cohomology equations for the zero average part. If we 
do so and substitute in the equation for W 2 , we are led to a linear equation for a when we 
impose that the averages of both sides in fl5.7p match. 

To do this computation explicitly, let us take the average of fl5.7p . which leads to solving 
the following equations for W 2 and a: 

0 = --SW2-W^-^{^a-A-Ma 
[X — 1) W 2 = —E 2 — A2(t , 

where Ba and Bb are such that (5a)° solves the equation \{Ba)^ — (i?a)° ° 
while solves the equation X{Bb)^ — o 7L = —(^ 2 )°. Let us write W 2 

for some unknowns A, B. Then, from the second of fl5.8l) we obtain 


(5.8) 

= 

= A + Ba 


A 


E 2 

A-1 ’ 


B 


A 2 

A-1 ’ 


which, substituted in the first of (I5.8p . gives a. 

In summary, we are led to the following algorithm (which is identical to Algorithm 33 of 
|CCd]L13c] L 


Algorithm 16. Given Ka : ^ Ad, /Xa G let X & C be the eonformal factor for the 

mapping f^^^. Perform the following computations: 

1 ) E^ f^^,,oKa-KaOT^ 

2) Q! •<— DKa 

3) N ^ [a^a]-^ 

4) M ^ [a I J-i o KaaN] 

5) 

6 ) E^(3E 

7) P^aN 

7 •(— a~^ J~^ o Kaa 

S^^iPo O lU J-' o KaP - X{N O o T^){N o T^) 

A^iMoT^)-^ 

8 ) {BaT solves A(i ?„)0 - o = -(^ 2 )°, 

{Bbf^lves X{Bbf - {Bbf o = -(^2)° 

9) Find 14^2, cr solving 

0 = -SW2-S{^-S{^a-¥,-Xa 
(A — 1)1T2 = —E 2 — A2a 

10 ) {W2f = (Ba)^ ^aW 

11 ) W 2 = {W 2 )° + IT 2 

12) (lTi)° solves (lTi)° - (lTi)° oT^ = -( ^ 2 )° - (.^ 1 )° - (Ai) V 

13) + MW 
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/4a ^ /4a + Cr . 

Remark 17. (/) It is important to note that AlaorithmITR involves only algebraic oper¬ 

ations, compositions of derivatives and solving cohomology eguations which work just 
as well when some of the objects involved are complex. Indeed, in |CCdlL13^ - and 
in good part of KAM theory - many functions are defined in complex extensions of 
the tori. 

(a) We note that, besides being the basis of the theoretical treatment in |CCdlL13^ . this 
is also a very practical algorithm, since each of the steps are obtained applying stan¬ 
dard algebraic manipulations, common in Celestial Mechanics. Notice that one step 
achieves guadratic convergence, but it is reguired a low storage and low number of 
operations. 

A consequence of Remark [T7I is the following result. 

Corollary 18. Consider that we give as input to the AlaorithmlTSi the members /C)/4e of a 
family indexed by e, with e ranging in a domain. Assume that the non-degeneracy assumption 
holds in the domain. 

If Ke, is an analytic (respectively continuous) family, then the result of the algorithm is 
also an analytic (respectively continuous) family. 

To obtain estimates for the iterative step described in Algorithm [161 we observe that the 
bounds for the correction {W,a) remain very similar to the estimates in |CCdlL13^ . The 
main difference with the procedure in |CCdlL13^ is that in Step 8, we use the estimates of 
Lemma [HI 

The estimates for the error in the step do not need any change from the treatment in 
[CCdlLblcj . since they just involve adding, subtracting, using the second order Taylor esti¬ 
mates and estimating the neglected term involving RW. 

After the estimates for the step are performed, we only need to check that the iteration 
can proceed and yields the desired results. This is nowadays quite standard and does not 
require any changes from the presentation in |CCdlL13c] to which we refer the reader for 
full details. 

We conclude by mentioning that the proof of Lemma [15] about the uniqueness of the 
solution does not require any change from Theorem 29 in jCCdlTdSc] . 


6. Prooe oe A) IN Theorem [12] 

In this Section, we prove part A of Theorem [12] which amounts to showing the existence 
of Lindstedt series to all orders. Let us start from the exact solution {Kq,p.q) as in (14.Ih : 
since we assumed that is symplectic, we have that 


■ 
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Let Mo = [DKq I J-1 O KqDKqN], N = {DK^DKo)-^ and let = <^ with ^ as in flOll : 
then, one obtains 

D/„,o o A'„(9) M„(9) = M„(9 + c^) ( “ ) , (6,1) 

Let Kj, inserting these power series expansions in the 

invariance equation fl3.ip . expanding the series in e and equating the coefficients of the same 
power of e, we obtain recursive relations dehning Kj and /ij, as described below. 

At the hrst order in e we obtain the equations: 

(L>/^g,o o Ko)Ki - Ki oT^ + o /Lo)/ni + o Kq = 0 , (6.2) 

while at the j-th order in e, 2 < j < N, we get: 

(-D//io,o ° Ko)Kj — Kj o o o Riii)fij = Fj{Ko ,..., Kj_i, po, fij-i) , (6.3) 

where Fj is an explicit polynomial in its arguments with coefficients depending on the deriva¬ 
tives of computed at /i^ = 0, e = 0 and composed with Kq. We note for future reference 
that 

1 F 

Fj{KQ ,..., Kj_i, po,..., /ij_i) = ° -^0 + ^j{Ko ,..., Kj_i, fiQ, ...,fij_i) , 

where Fj does not depend on derivatives of order j of Equation fl6.2p is of the same 

kind of (I6.3p . since it suffices to dehne Fi{Ko, hq) = o Kq. 

To solve fl6.3p (equivalently fl6.2p L we write Kj{6) = Mq{ 6)Wj{6) for a suitable function 
Wj = Wj(9), so that fl6.3p becomes 

o Kq) MqW j — {MqW j) o -|- {D^f^^ Q o iLo)hi ~ ■••i hoi ■••i hj-i) • 

Using fl6.1ll we obtain: 

( 'o ^Id ' ) ° ° 

= [MqoT^) ^ Fj{KQ,Kj^i, fiQ,^j_i) . (6.4) 

Writing fl6.4l) in components, again we recall that this means taking the hrst d rows and the 
last d rows, say Wj = {Wji\Wj 2 ), we get the following equations: 

W,2-Wj20T^ + [{MQ0T^)-\D^f^,^Q0KQ)]2^i, = [{MqoT^)-^F,]2 
o + SQWj 2 + [(Mo o o Kq)], fi, = [{Mq o T^)-^F,], ,(6.5) 

where [-ji, [•]2 denote the hrst and second component. Under the non-degeneracy assumption 
HND and provided to G Prf(i^, r), equations fl6.5p can be solved to determine hU,i, hU,' 2 , jOj, 
according to the following procedure. Let 

Ej = (Mo o T^)-^F, , Ao = (Mo o T^)-\DJ,,,q o Kq) ; 

then, fl6.5p becomes: 

Wj2 — Wj2 oT^ + A2oP-j = Ej2 

Wji — Wji oT^ + AioPj = Eji — S'ohU,'2 . 


( 6 . 6 ) 
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Taking the average of the hrst equation in (16.6p . we obtain 

^20 hj = Ej2 , (6-7) 

which determines fij. Taking the average of the second equation, we obtain 

= Eji — {SqW j2) ■ ( 6 - 8 ) 

Let Wj 2 = Wj 2 + (hTj 2 )°; then, we have 

{SoW^2)=^oW~2 + So{Wj2f . (6.9) 

Using that Wj 2 is an affine function of yUj, we write {Wj 2 )^ = (-Bao)° + where (i?ao)°, 

(i?bo)° the solutions of the equations: 

{Baof-{BaofoT^ = {E,2f 

{B,of-{Bu,foT^ = -{A2of . ( 6 . 10 ) 

From fl6.8l) . (16.91) and (16.101) we obtain: 

+^oW~ 2 + SoiBto)^ /i, = -So{Bao)^ + E~i . (6.11) 

The equations (16.7p and (16.1 ip in the unknowns Wj 2 and /ij can be written as 

Sq Aiq + 5'o( i?feo)° A / Wj2 \ ^ —5'o(i ?ao)° + Eji 

0 A 20 y V h-i / Y Ej2 

which can be solved to obtain Wj 2 , fij, provided HND is satished. The proof is completed, 
once we solve the equations (16.61) for the non-average parts of Wi and W 2 . The equations 
(16.61) are cohomo logical equations of the form (12.8p with A = 1. 

Assuming that we have determined the functions Kj and the terms fij for 1 < j < iV, 
we obtain the hnite sums kI~^\ which solve the invariance equation within an error 

given in (14.2p . 



7. Prooe oe B) in Theorem O 

We start by considering an approximate solution, as provided by part A of Theorem [I2l 
namely a solution (A:|-^\/i'-^'), which can be expanded in formal power series in £ as 

N N 

j=0 j=0 

and which satishes the bound (14.2p . 

Let A > 0 and let £0 ^ ^ro(^) with the set Qro as in (12.6p . where the cohomological 
equations can be solved. Assume that e belongs to a sufficiently small ball B centered in Eq, 
for example 

k-^ol < 10“® dist(eo,dQro(^)) ■ (7-1) 

By the choice of e in this subset of ^ro(A), all the assumptions stated in Theorem [14] are 
satished. In particular, the Diophantine condition HI is required also in Theorem [T2j The 
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approximate solution in H2 is provided by the choice {Ka, f^a) = the associ¬ 

ated error term can be assumed to be sufficiently small as in H5 thanks to the inequality 
fl4.2p and the assumption flT.ip . provided tq is sufficiently small. We recall that due to HA 
we have A(£) — 1 = -f 0(|e|“+^) for some a G « G C, so that the right hand side of 
fl4.2p can be estimated by a power of |A(e) — 1| as 

° Ah"' - Ap' °rjlp. < C'^\x- 1|™ 

for some constant C'^ > 0. By taking A close to 1, we obtain the bound in H5 on the error 
term. 

Using the fact that the determinant is a continuous function, the non-degeneracy assump¬ 
tion H3 is implied by HND, provided HA is satished and tq is sufficiently small. 

In conclusion, all assumptions required in Theorem [H] are satished and Theorem fTTI allows 
one to state that there exists an exact solution of the invariance equation fl4.3p . 

The conditions of Theorem [14] are verihed uniformly and therefore the sequence of ap¬ 
proximate solutions constructed in the proof of Theorem [T4| converges uniformly to the true 
solution (i^£,/ie) satisfying (14.31) . Such families of functions will be analytic in e in the 
interior of Qro and continuous in all of Qro as a consequence of Corollary [181 

Due to the construction of the exact solution, the inequalities fl4.4l) are satished. 

We conclude by mentioning that the error is as small as we want for e small and N large 
enough, and that the bounds on the constants depend on the dehnition of the set fl7.ip . 
As remarked before, we can ensure that the solution satishes the normalization (13.41) : such 
normalized solutions are locally unique. 

8. Further geometric properties of the sets Q , A 

In this Section we formulate two geometric properties of the set A dehned in (12.7p . From 
this, one can obtain properties of the set Q dehned in fl2.5p . since it is obtained from A by a 
conformal transformation. 

As a matter of fact, we will also consider the set 

T(A; (V, r) = {A G C : z/(A; cn, t) < A} , 

which is easier to study. Clearly we have that in a neighborhood of the unit circle 

A{A-,u,t,N) D T(2.1“^^’''^M, (n, r) . 

This inclusion is a not very accurate estimate near £ = 1, since we are ignoring the factor 
|A — 1|^+^. On the other hand, for properties over the whole unit circle, this is a good 
estimate. 

We observe that, since the number u is the supremum of several quantities, the sublevel 
sets are obtained by removing the sets where one of the inequalities required in the dehnition 
of T fails. We observe that the places where one of the inequalities fails can be bounded 
from above by a ball and we can obtain a lower bound for the set T by removing from the 
plane balls which enclose the region where one of the inequalities fails. These inequalities. 
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and hence the excluded balls, are indexed by /c G Z'^\{0}. The radii of the excluded balls 
decrease as |/c| grows. 

Of course, the set T contains open balls outside the unit circle. Hence the interesting 
question consists in studying the properties of density of T on the unit circle. 

In this Section we will establish two geometric properties of the set A: the hrst one 
states that one is a point of density for the set A, both as subset of the complex and also 
restricted to the unit circle (see Proposition [2T]) : this result is independent of A. The second 
property (see Proposition [20]) shows that the set of points that are tangentially accessible 
(see Dehnition [19] and compare with |Car54b] . |Car54a] . |Car52] i is also of large measure 
near one. Both results are proved by the standard argument in Diophantine approximation 
theory by estimating the measure of the excluded balls. 

Later on, in Section [10.41 we will show that the set Q can be improved to another set Q, 
which is tangentially accessible in more points (compare with Proposition [2^ . 

Definition 19. Let C be a complex domain. We say that a point Aq E C is tangentially 
accessible in C, when there exists a unit complex number u (denote by it its complex conjugate) 
and there exist d > 0, 'j > 0, m > 2, such that 

P(Ao, m; 7, (5) = {Ao + + sh : |f|, |s| < 5, s > 7|f|™'} C C . 

We say that a point is tangentially accessible from both sides in C, when 

P^(Ao, m; 7 , (5) = {Ao + tu + su : |f|, |s| < 5, |s| > 71^1™} C C . 

When we need to me more precise, we can talk about an m-tangentially accessible point. 

Note that, clearly if a point A G C is tangentially accessible for a set T and T C A, then 
A is tangentially accessible for A. 

The property of being tangentially accessible is, of course, only relevant for the points 
in the boundary. It means that we can get regions bounded by parabolas tangent to the 
point inside the domain. If the boundary is given by a differentiable curve, all the points 
are tangentially accessible. Also, when we transform a domain by a differentiable mapping, 
all the tangentially accessible points for the original domain get mapped into tangentially 
accessible points for the image. 

The fact that a point is tangentially accessible has important consequences. For example, 
in [SokSOl IHar49] it is shown that the asymptotic expansions based at that point can be 
Borel summed and determine uniquely the function in the sector; the book |Car54aj contains 
several other properties, which are a consequence of accessibility. 

In our case, the points in the boundary of A are tangentially accessible from both sides. 
We also note that the asymptotic expansions constructed in Section [6] are dehned on both 
sides of the domain. It then follows that we can continue the functions in a unique way 
across these points. 
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Proposition 20. Consider u G a) and assume that for some m > 2 the following 

inequality holds: 

a > md . (8.1) 

Then, almost all points in the unit circle are m-tangentially accessible for UaT{A]u, a). 
In particular, by taking A sufficiently small, we can get a set of m-tangentially accessible 
points, whose complement has measure as small as desired. 

Proof. Remember that an upper bound for the set T is obtained by removing balls centered 
at of radius approximately equal to A~^\k\~'^. 

We observe that for all the points in the unit circle that are at a distance bigger than 
CA~^\k\~^Y^'^ from for some positive constant C, we can hnd a domain of the 

form r(a,m), which does not touch the excluded ball centered at 

Therefore, for all points in a subset of the unit circle whose complement has measure less 
than the condition imposed by the ball corresponding to k is not an 

impediment for being m-tangentially accessible from both sides. 

It then follows that the set of points that satisfy the conditions imposed for all k to be 
m-tangentially accessible from both sides has a complement whose measure is less than 

fceZ'^\{o} teN 

for some constant C"; the factor 2 at the l.h.s. takes into account that we consider both sides 
of tangential accessibility. We see that under the condition (18.Ih . the sum above is hnite 
and, by choosing the constant 7 large enough, we can obtain that the complement of the 
m-tangentially accessible points in A(y4; uj) has measure as small as desired. □ 

Proposition 21. Assume that uj G t) with 2 t > d and that N > 1. The point X = 1 is 
a point of density for the set A G C (with the two-dimensional Lebesgue measure). If r > d, 
then 1 is a point of density for A fl (with the one-dimensional Lebesgue measure). 

Proof. This is a standard excluded measure argument. Fix p > 0 sufficiently small. In the 
following, we will not specify the constants that we generically denote as C. Let Cp = {A G 
C : p < |A — 1| < 2 p}. 

The set A^ U Cp is the union of the sets 

7^fc,p = {A G C : - A| < A-^\k\-^\X - l|^+\ p < |A - 1| < 2p} . 

Note that, in particular, if there is a point A satisfying the two conditions dehning IZk^p, we 
have 

- 1| < 2^+^ p^+^ + 2p < 3p , 

where the last inequality holds for p small enough. Since uj is Diophantine, this implies that 

\k\ > Cp--r. 

Now, for a hxed k we see that the set of points A for which A| < 2^~^^ A~^p^~^^\k\~'^ 

is a circle with area smaller than CA~‘^p'^^^~^^'>\k\~‘^'^. 
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Hence the total area excluded is less than: 


k&d.\^0},\k\>Cp-r 


^j^-2p2{N+l)\^-2r 






We observe that, under the hypothesis 2t > d the measure of the regions excluded in an 
annulus of inner radius p and outer radius 2p is less than 

^^-2p2{N+l)+^{2r-d) _ 


Hence, the excluded measure in the ball of radius p around the origin is bounded by a power 
greater than 2. Indeed the power is arbitrarily large if N is large enough. 

The argument for the measure excluded in the unit circle is extremely similar. The only 
thing we have to change is to estimate the length of the excluded intervals rather than the 
area of the balls. We obtain an estimate for the excluded length less than 


□ 


9. A CONJECTURE ON THE OPTIMALITY OE THE RESULTS 

The domain of analyticity of the embedding and the parameter can be optimized 
in a variety of different ways; for example, the constants in the cohomology equations can 
be sharpened, the order of the expansion can be taken to optimize the results, by exercising 
more care in the method presented here or, perhaps, developing a new method, etc. Indeed in 
Section no.41 we will present an improved domain of analyticity. Nevertheless, in this Section 
we want to argue that the shape of the domains established in Section 12.2.11 is essentially 
optimal in the following sense. 

Conjecture 22. Consider a generic family of mappings f^^^, satisfying the hypotheses of 
TheoremlTB. Consider the solution {K^, pf) produced in TheoremlTB and the maximal domain 
where such solution is defined. Then, there exists a seguence of balls Bk in the complex plane, 
which are not included in the maximal domain of analyticity. 

The balls Bk with centers in Sk, such that X{ek) = correspond to the balls excluded 

in the definition of Q (see fl2.5p ). 

We note that if Conjecture [22] were true, it would have important consequences for the 
analytic properties of the asymptotic expansions. 

The conjectured analyticity domains do not contain sectors centered at the origin with 
aperture bigger than n/a. When a > 1, we cannot use the Phragmen-Lindelof principle 
f |PL08j . |SZ65j . [Har49 ] ) to obtain that the Taylor asymptotic expansion determines the 
function. Indeed, one could get non-trivial analytic functions with zero asymptotic expan¬ 
sion. It is not clear whether any method of summation based only on the expansion produces 
the correct result solving the functional equation. 
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The main reason for stating Conjecture [ 22 ] is an argument which goes by contradiction, 
already used in |CCCdlL15] to which we refer for more details. We recall that the drift /i, 
together with the embedding, is an unknown of the problem. Then, we extend the one- 
parameter family e, depending on the parameter £, to a family depending also on a 
parameter 7 , such that 

Then, we proceed as follows: 

{i) we assume that is analytic in the parameters £, 7 ; then, under suitable resonance 
conditions on the conformally symplectic factor (see Proposition [23] below), that 
ensure that A(efc) = for an inhnite sequence k e Z'^\{0}, we claim that 

there is no family of functions analytic in 7 for 7 small near 
{ii) by recalling a result stated in |CCCdlLl^ . we show that if we have solutions in 
the space of analytic functions for all in a neighborhood, there has to be a family 
analytic in 7 for 7 small; 

{Hi) the consequence of the above two statements is that in every neighborhood, there 
has to be a family which is analytic in the neighborhood of the resonant points 
where X{ek) = 

The result in {i) is given by the following Proposition. 

Proposition 23. Let the family of mappings f^^^ satisfy the hypotheses of Theorem [TM 
Assume that there exists an analytic solution {K^,/j,^), solving the invariance eguation (13.Ih . 
Let fij.,e,'y be a generic family, which is analytic in e, 7 , with conformally symplectic factor 
A = X{e,'y) and such that 

fli,e,0 ffi,e ■ 

Then, we can find an infinite seguence Sk, k G Z'^\{0}, with X{ek) = such that there 

is no formal expansion in 7 for the solution of fl3.ip associated to 

The proof of Proposition [23] is given in Section 19.11 

Notice that the proof of Proposition [23] shows that near e = Sk, the effects of a change 
on {K^,fi^) are much larger than the change on Hence, it seems likely that one can 

destroy the solution {K^, fie) without altering too much the family 

We interpret Proposition [23] as meaning that the solutions for families which have large 
domains of analyticity are unstable and can be easily destroyed. Of course, the sense in 
which we prove instability is somewhat weaker than what is needed to reach the conclusions of 
Conjecture [22] rigorously, but it goes in the right direction and this motivates the formulation 
of the results as a conjecture. 

We note that the above argument suggests that the set of functions such that 

has a torus that is analytic in a neighborhood of Sk, is a set of inhnite co-dimension (in 
particular nowhere dense). Once this result were established, the set on which inhnitely 
many resonances are destroyed is residual, in the sense of Baire category. 
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Note that our argument in Proposition [23] is somewhat similar to the arguments in [Poi87] 
about the lack of uniform integrability; indeed, what we conjectured is an analogue of the 
lack of integrability for generic systems, based on the lack of uniform integrability. 

9.1. Proof of Proposition l23l The proof of Proposition 1231 is very similar to the results we 
obtained on Lindstedt series in Section [ 6 ] If there existed an approximate solution 
satisfying 

O Ke- KeOT^ = E 

with an error E, then we can hnd a matrix M such that 

o K { 0 ) M { 0 ) = M {0 + !.,)(' “ ffj )+ R , 

where M and S have been defined, respectively, in fl3.2l) and fl5.3p with ^ replaced by 
//ie,£, 7 ) R is the error in the automatic reducibility as in fl5.5l) . By the procedure of 
automatic reducibility, we see that after making expansions in 7 , we have to solve the 
following two equations: 

Wi-WioT^ = -sw2-A-M(y 
AIT 2 — W 2 o Ttj = —E 2 — ^2(7 , (9-1) 

where W = (hPi, 1 ^ 2 ) and a are the corrections to the step, E = {MoT^)~^E, A = [Ai|y 42 ] = 

We must require that the right hand sides of (19.11) have zero average. 
This can be obtained by properly choosing a, provided that the non-degeneracy condition 
(15.21) is satisfied. 

We note that precisely at the points Sk where X{ek) = k G then the 

second equation in (19.1!) has obstructions to solutions. We also note that, because of the 
form pointed out in (19.11) . the solution of the order j equations of a formal expansion in 7 
exists only if ■£^fij,^,e,'y\'y=Q satisfies a condition. Of course the conditions for different j are 
independent, since they affect different coefficients. Therefore, the existence of expansions 
to order j requires that the mapping belongs to a co-dimension j submanifold of maps. □ 

As indicated in the preliminaries, we note that for the values £ ~ we have that the 
changes on fj,^ induced by perturbations of are much larger than the perturbations 

of ffie,£n itself, this makes it plausible that one can introduce changes in ^ which destroy 
the analytic without destroying (see |Sie54j for similar arguments). 

10. Automatic reducibility and Lindstedt series 

In this Section we present several results about Lindstedt series. First, we show that 
we can obtain Lindstedt series expansions around any point £0 (see Section riO.ll) . Some 
relations with the theory of monogenic functions are presented in Section IIP.31 Then, by 
lifting the automatic reducibility to families, we show how to get quadratically convergent 
algorithms for the Lindstedt series; as a byproduct of this result, we will obtain an alternative 
proof of Part B of Theorem [12] (see Section [10.2p . Finally, we will show how the results of 
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Section 110.11 lead to define an improved domain of analyticity (see Section llO.dp and we 
conclude by establishing the Whitney differentiability of K^, on Q (see Section flO.Sp . 


10.1. Lindstedt series from any analytic torus. In this Section we show that if for some 
Eq there are solving the invariance equation at £o) and which satisfy some mild non¬ 

degeneracy assumptions (which are implied by |eo| ^ 1), we can hnd a formal power series 
in e — ^0 that solves the invariance equation in the sense of formal power series expansions. 

Of course, for points in the interior of the analyticity domain, the existence of expansions is 
obvious. The interesting case is that the same result holds for some points in the boundary of 
the analyticity domain. As we will see, this matches very well with the theory of monogenic 
functions (see Section FlO.hp . It is also important to notice that even in the interior of the 
domain of analyticity the present method gives very good estimates of the formal power 
series, much better than what can be obtained from just Cauchy estimates, see Remark [27l 


Proposition 24. Let u be a Diophantine vector in the sense of Definition 0 let f^^i. he a 
family of conformally symplectic systems as before. Assume that for some Eq we can find 
G Ap, /ieo G such that = A(£o)f^ o,nd 

/^eo.eo ° Keo = O . (10.1) 


Assume furthermore that A(£o) is Diophantine with respect to uj in the sense of Definition^ 
that iC£g(T(() is well inside the domain of definition of and that satisfies the non¬ 

degeneracy condition H3 of Theoremll^ 

Then, for any 0 < p' < p, there is a formal Lindstedt power series solution kI°°\ 


74=0 

CO 

74=0 


( 10 . 2 ) 


with coefficients Kn G Ap', fin G that satisfy the invariance equation in the sense of 
formal power series, namely 




O Ad<A^] _ o T„ 


^ r< \c c- I'^+i 

A — ^oI ) 

p' 


(10.3) 


where = Y.n^hn(e - eg)". 

The coefficients may be chosen so that the normalization (13.4p . with M the corresponding 
for is satisfied in the sense of formal power series. The series that satisfy fllO.ip and 
(13.41) in the sense of power series are unique. 


Proof. Introduce the matrix M corresponding to Keo (13.21) . for which we have 


CT.„,..oA-..Af(«) = M(« + ij) 


Id S(ff) \ 

0 A(e„) Id J ■ 


(10.4) 


Since the invariance equation has zero error, so does the reducibility equation fll0.4p . 






















34 


R. CALLEJA, A. CELLETTI, AND R. DE LA LEAVE 


If we substitute fll0.2j) in the invariance equation and equate the coefficients of order 
{e — on both sides, we obtain for n > 1: 


£’//.eo,£o ° KeoKn + ° ^eo/^n + Rn = Kn O , (10.5) 

where Rn is a polynomial expression in ii'i, ... K^-i, /ii, ... /in-i with coefficients which 
are derivatives of / evaluated at 

We think of fllO.51) as a recursion that allows one to compute Km fin, once Ki,... Kn-i, 
yUi,... fin -1 have been computed. 

We also substitute (I10.2p in the normalization condition (13.dh so that we can obtain expan¬ 
sions of functions that satisfy the normalization condition with respect to We obtain 
that the coefficient of order [e — So)^ of the normalization is just 


ijd 


M-\e)Kn{e) 


de = Q 


( 10 . 6 ) 


Of course, the equation (110.51) is a particular case of the equations appearing when we 
studied the Newton’s equation for approximately invariant solutions in Section [50] fit suffices 
to take E = Rn). Many more details can be read in [CCdlLlO^ . We note that using 
the change of variables Kn = MW reduces the equation fll0.5p to constant coefficients 
cohomology equations. We emphasize that the series whose coefficients satisfy (llO.Op are 
unique. □ 


Two important corollaries of Proposition [2l] are given below. Corollary [25| shows that one 
can construct asymptotic series of the form (I10.2p . which in principle do not converge, such 
that the truncated series to order N satishes the invariance equation up to a small error. 
This result allows to start an iterative process from the approximate solution, leading to an 
exact solution which is locally unique. 

Corollary 1261 shows that the derivatives of fi^ are obtained by hnding the expansions; 
their estimates involve only the loss of domain in the coordinates, while they are uniform in 


Corollary 25. Within the assumptions of Proposition let Eq he such that we can find 
G Ap, fisQ G C'’* satisfying (llO.ip . For any 0 < fi < p, let kI°°\ pt^^he such that (110.21) 
holds with coefficients Kn G Ap, Pn G and that the truncated series p^^^^ satisfy 

(llO.Op and the normalization condition. 

Then, for all e sufficiently close to Eq, one has 


N 


J2Kn{e-Eo)^-K, 


n=0 


P' 


N 


'^Pnie - SoY 


Pe 


n=0 


< Ci,[\e — Eq 


< Cn\s — Eq 


V+1 


\N+1 


(10.7) 


for some positive constant Cn depending on p'. 
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Proof. From the construction of the series we obtain that the polynomials kI~^\ satisfy 
the invariance equation up to a small error as in fllO.Sp . Hence, all assumptions H1-H5 of 
Theorem [14] are satisfied. Applying Theorem [TH with Ka = kI-^\ /i^ = we obtain an 

exact solution K^, pg, which satisfies the bounds fll0.7p . Using the local uniqueness of the 
normalized solutions in Lemma 1151 we obtain that the solution produced is K^, □ 


Corollary 26. Assume that for e E Q we have 


Ke\\p < B, 
\Pe\ < B 


for some B > t). Then, we have for 0 < p' < p: 

< Cj{p- 




d y 

Te) 


< Cj{p - B , 


( 10 . 8 ) 


where the derivative is understood in the regular sense for e in the interior of Qro ond as the 
(unique!) term determined from the Lindstedt expansion; the constant Cj depends also on 
the Diophantine constant. 


Proof. At each step, the derivatives in (110.8p are obtained by solving two cohomology equa¬ 
tions: one of them is a regular cohomology equation not involving the conformal factor, while 
the other depends on A (see Proposition |24| and compare, e.g., with (19.ip ). At each step one 
has a constant loss of domain (p — p'). The solution of each equation gives a factor which 
contains a power of the loss of domain (p —as stated in Lemma IH1 The constant Cj 
depends also on the quantities appearing in the definition of the set Q in (12.5p and precisely 
on the product z/(A;a;,r). □ 


Notice that we have stated that the bounds of | (^)'^ Pe\ depend on p — p', because this is 
what comes directly from the proof. Of course, ps does not depend on p, so that we could 
simply get 


dy 

Te) 




Remark 27. It is interesting to compare Corollary\2^with the results obtained by applying 
the Cauchy bounds by drawing just a small circle. Of course, Cauchy bounds produce bounds 
on the derivatives in the same space Ap and the bound involves the distance to the boundary. 
The bounds in Corollary are uniform in the distance to the boundary, but are only true 
in a slightly weaker space. 

The main reason for the difference in the results is that we use not only function theoretic 
properties, but we use that the functions we are studying satisfy a functional equation. The 
use of the functional equation to obtain the derivatives provides much better estimates than 
those available just using the function theoretic properties, albeit in a different space. 
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Of course, there are many functions that are defined in the same domains, but which, not 
satisfying any functional equation, have much worse properties. The systematic use of the 
functional equation and the fact that this functional equation leads to automatic reducibility 
is one of the big advantages of the present method with respect to using only function theoretic 
methods. 

10.2. A quadratically convergent algorithm for Lindstedt series and an alternative 
proof of Part B of Theorem 1121 An alternative way to compute the step by step series as 
in part A of Theorem [12] is to lift Algorithm [16] to the family of maps depending also on the 
parameter e. Furthermore, the lift of the algorithm leads to a proof of Part B of Theorem [12] 
based on an implementation of Newton’s method after introducing a lift to functions of e 
and performing a Newton’s step as a function of 9 and e. 

In the following, we show that lifting the algorithm to compute families of functions 
depending on the parameter £ is a very natural procedure. We consider the family of 
conformally symplectic maps. Assume that for some e such that \e — e^ol is sufficiently small, 
we start with an approximate solution 

N N 

Ki^N] ^ ^ ^ ^[<iV] ^ ^ ^ 

n=0 n=0 

which satishes the invariance equation up to order N + 1, i.e. 

/ [<.] ^ O O (10.10) 

H-e 

with the following bound on the error for some 0 < p' < p; 

l|£P'IU<C|e-£„r+' , 

Let us also dehne: 

N 

n=0 

Algebraic manipulations similar to those performed in the proof of Theorem [14] lead to the 
expression: 

° + £,) ^ “ ) + me ) (10,11) 

for a suitable error function R = R{9), where (again) we dehne the functions mI~^\ Ne~^'^ 
as 

M,{e) = [DKl-^\e) I j-^ o Kl-^\e)DKl-^\e)Nl^^\e)] 

with 

and as in fl5.3p . We look for a new approximate solution, given by the updating of 

the approximate values and with corrections (equivalently Wf) and i.e.. 
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K'^ = + Ae = /ig = + (Je. Proceeding in the same way as in 

Section [5l we look for the corrections (hhgjCTg), that satisfy, 

('o ap'i Id ) ° " “ (Mpior„)-TP' 

- (Mjs"! o o irp^V.. 

which in components provides the following cohomological equations similar to fl5.7p : 

W^-WioT^ = - El - IiCT 

XfN]w^_W^oT^ = -E^-Aa, ( 10 . 12 ) 

where E = {E,, E^) = oT^)-^ E'f^\ A = [I 1 II 2 ] = ^ojppi). 

M'S ;C 

The functions E, A, Wi, and IT 2 depend on the parameter e, but we omit it in the notation 
to avoid cluttering. 

We notice that the cohomological equations in fll0.12p can be solved under the analogous 
of the non-degeneracy condition fl5.2p . Then, the solution of fll0.12p provides the desired 
corrections (H4, a^). We emphasize that the objects involved in the solution of the cohomo¬ 
logical equations in fll0.12p are functions of e and that the estimates used for the iterative 
step in Algorithm [16] are uniform in e. Furthermore, if we use the supremum norm the same 
estimates for the iterative step hold for the convergence in the families (see |dlLQ00] j. Of 
course, with this sup-norm the convergence is uniform for the whole e dependent family. 
Since the uniform limit of analytic functions is analytic, the solution must depend analyti¬ 
cally on e. 

The iterative method also suggests a practical algorithm for the computation of the per¬ 
turbative expansion in e. If we start the iterative procedure from an approximate solution 
having terms up to {e — eo) to the A-th power, the quadratic convergence of the Newton’s 
method implies that after every step the number of terms in the expansion doubles. One can 
use methods in automatic differentiation (see [Harlll IBK78[ IBis09] l to implement the oper¬ 
ators involved in the iterative step, so at every iteration the new step produces N new terms 
and it also has the advantage that it corrects the previously computed lower order terms. 
Furthermore, on top of being quadratically convergent the Newton’s method is numerically 
stable. 

Moreover, if A satishes the Diophantine assumptions, within the regions in e where we 
have good estimates one can hnd a sequence of analytic functions converging uniformly in 
the domain Q^o- In the set Qro we obtain the convergence of the sequence of approximate 
solutions, obtained applying the iterative steps, to the exact solution with bounds analogous 
to those in fl4.4p . Hence, the fast algorithm described before provides also a different proof 
of Theorem [121 

10.3. Relations of the previous results with the theory of monogenic functions. 

We recall that a function k dehned on a set "H C C is monogenic, when for all Eq in Ti, the 
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limit 


lim 


k{e) - k{eo) 


(10.13) 


e-So 

exists. 

When l-i is an open set, this is the nowadays current dehnition of differentiable functions, 
but for more complicated sets, since the limit in fllO.ldp is taken only when £ G "H, the notion 
may be very different. The book |Borl7j contains an account of the history of monogenic 
functions and their relation with the nowadays standard notion of differentiability. A problem 
of great interest is to study the relation between monogenic functions in some sets and quasi- 
analytic classes (that is, classes of functions such that the function is determined by its values 
in a set). The modern standard result establishing that analytic functions in an open set 
are determined by their values on a set with an accumulation point is a predecessor of those 
results in which the structure of the set determines the function. It is also known that 
the monogenic properties in domains with enough aperture lead to consequences for the 
properties of asymptotic expansions f |Win93] b 

The result stated in Proposition [2l] implies in the classical language that the functions 
are monogenic at many points (namely, the points for which X{e) is Diophantine with 
respect to u) in the set Q. 

The relations of KAM theory with the theory of monogenic functions has been studied 
in |Kol571 IArn61( IHer85( IMSn3( IMSlH ICMS14] . These papers investigate the monogenic 
properties of the solutions, when the frequency u goes into the complex. When u has 


complex coefficients, then e 


2'Kik-U) 


may be bounded away from 1. Hence, the cohomology 


equation does not present any small divisors and, indeed, it may be a compact operator. 

The approach of the above papers to study the monogenic properties is very different 
from ours. They study in great detail the monogenic properties of the solutions of linearized 
equation and use this as the basis of an iterative method. In our case, we just use the 
formal power series and the a-posteriori format of Theorem [HJ We take advantage of the 
non-perturbative properties of Theorem [14] and the fact that the solutions are characterized 
by solving a functional equation. 


10.4. An improved domain of analyticity. In this Section we will show that the domain 
Q can be slightly improved. This does not contradict Conjecture [22l since the extensions 
constructed here are precisely in the places away from the resonances which Conjecture [22] 
claims are essential. 

The following Proposition [28] is the main result of this Section. We show that one can hnd 
a larger analyticity domain for iC, p in such a way that the new domain is m-tangentially 
accessible at any Diophantine point (the analyticity of K should be understood in a slightly 
smaller domain). 

The proof of Proposition [2H] is based on Proposition [24] that established that we can make 
a Lindstedt expansion around any Eq sufficiently small provided that A(£o) is w-Diophantine. 
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Proposition 28. Under the assumptions of TheoremlTM and Proposition^^ let the functions 
K : Qro ^ U ■ Gro ^ solvc fl3.ip . 

For any 0 < p' < p, we can find an analyticity domain G, such that any point eo sufficiently 
small is m-tangentially accessible, m E N, from both sides in G, provided Sq is such that 
|A(£o)| = 1 and A(£o) is Diophantine with respect to oj. 

The functions K, p are analytic from G to Ap, TA, respectively. 

Proof. To prove Proposition |2H1 we observe that, given any point £o as in the assumptions, 
using Proposition |24] on the existence of Lindstedt power series in £ — Eq, for any N we 
obtain a polynomial that solves the invariance equation up to an error measured in the p- 
norm which is smaller than Cjs[\e — Precisely, we construct the polynomials = 

Yln=0 - SoT, = J2n=0 “ ^o)”, SUch that 

11/ Kiv, ^ o iPpi - iPpi o T4p < CnIs - £ 01 "^+' 

P'S 

for some positive constant Cn- 

We also observe that, according to Remark [3, when we move out of the unit circle, the 
number A becomes cn-Diophantine with a Diophantine constant which is bounded from below 
by the distance to the circle. Precisely, if |A(e)| ^ 1, then we obtain 

z/(A(£);a;,r) < [l - |A(e)||"^ . 

If we fix 0 < p' < p and we denote 6 = p — p', we can apply Theorem [TTl provided that the 
error in H5 is sufficiently small, namely 

CnIs - \ l - \X{e)\\^ < Cu{uj-,Tf . (10.14) 

We now observe that, because A(e) is a conformal mapping which is at most tangent to 1 
of order a due to HA, and N is arbitrarily large, we can find a domain of e m-tangentially 
accessible (m arbitrarily large) where (llO.ldh is satisfied; then, we can find an exact solu¬ 
tion of the invariance equation fl3.ip for any e E G- Again we note that this solution can 
be obtained to be normalized with respect to (13. 4 p corresponding to K^. Then, the local 
uniqueness of normalized solutions tells that this should agree with the original solution. □ 

It is interesting to compare Proposition [28] to Proposition [201 Of course, they refer to 
(possibly) different sets: Proposition [20] is based just on a measure theoretic argument and 
only asserts the existence of points occupying a set of large measure, whereas Proposition [28] 
gives a concrete criterion for the points for which it applies. Of course, given the easy measure 
estimates of Diophantine numbers, there must be points for which both propositions apply. 

We note that the assumption that Eq is sufficiently small can be substituted with the 
requirement that there is a solution satisfying the invariance equation. We also note that 
the assumption |A(£o)| = 1 is, of course, not needed. When |A(eo)| 7^ 1, h is easy to shovj^ 
that we can get a solution defined in a ball around ^o- 

^Just note that for all e close to eo, we have that |A(e)| 1 and, in particular, A(e) is Diophantine. For 

any e close enough to eo we can take as initial point of the iterative step Kq, pq. 
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Remark 29. Even if Proposition\2B involves a loss of domain in the 6 variable, this is not a 
problem, because the regularity in the parameterization variable can be bootstrapped. Indeed, 
in jCCdlLLScj . there is a result establishing that all finitely differentiable solutions of the 
invariance eguations are analytic. 

10.5. Whitney differentiability of the parameterization Kg, and the drift on 

In this Section, we establish that the functions it'e, are Whitney differentiable in 
The argument we present here takes advantage of the fact that Ke,pie solve a functional 
equation. This allows us to obtain expressions of the formal derivatives with respect to 
parameters, since these formal derivatives satisfy functional equations which we can solve 
using the automatic reducibility. Using the a-posteriori format of Theorem [TU we obtain 
that the formal expansions lead to a small remainder. The rest of the argument uses that 
the domain Qr^ is a compensated domain, that means that the points can be joined by a 
path in the interior whose length is comparable to the distance between the points. There 
are many other similar problems in KAM theory and it is possible that similar strategies 
may work in them. 

Let us start by recalling the classical dehnition of Whitney differentiability f |Whi36tlSte70 
IFefn9[ IGraMj ). We note that this is a very real variable dehnition. In our case, when we are 
considering functions of a complex variable, we have, of course that C = 


Definition 30. Let T> C M"', n > 1, be a closed set. For i G Z_|_, t < p < t -L \, let H be 
a Banach space. We say that the function h : V ^ H is -Whitney differentiable (which 
we denote by h E C^f^fD)), when for every x E V we can find a collection of continuous 
functions indexed by k E N", say {h^x'^}o<\k\<i, taking values in X, such that 

||hi^)||<M, 0<k<i, xeV (10.15) 


and for all x,x -\- A E T> 


h 


(k) 

x+A 


E A' 

\j+k\<£-^ 


< M . 


(10.16) 


Then, the function h is called -Whitney differentiable in V with Whitney derivatives 
{h^x^}\k\<£- 


If P = M”, one can see that the case /c = 0 of fll0.16p is just the Taylor expansion of h, if h 
was differentiable. The cases /c > 1 of fllO.lOp are the Taylor expansions for the derivatives. 
It is known that when T> = M”, one needs to verify only the k = 0 case of fll0.16p . Then, 
the celebrated converse Taylor theorem f |AR67[ lNel69] j shows that the function h is indeed 
p-differentiable and that therefore, the j-th derivative is {p — j)-differentiable. 

When P is a small subset of M”, the Whitney derivatives do not need to be unique (think 
for example of the case when D is a hyperplane; since fllO.lGp only involves what happens 
in the plane, we see that we can put anything we want in the transverse directions). Hence, 
in general, the Whitney derivatives are not unique and one can see that the case /c = 0 in 
fllO.lGp does not imply the cases for higher k. 
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The remarkable result of |Whi36] is that given a function h which is C^-Whitney differ¬ 
entiable in a closed set "D, there exists a C^-function (in the standard sense) in M” which 
agrees with h on V. The original result is presented only for real valued functions, but it is 
easy to check the proof when the functions take values in a Banach space, since the proof 
follows from the proof in M” line by line (in contrast, the case when V is contained in a 
Banach space seems to be an open problem). 


Remark 31. There are some well known extensions, which we will not discuss here, but 
which we mention. In |Ste70] . one can find results showing that the extension can be obtained 
to be linear in h. Actually, one can construct a linear extension operator Ei which works 
for all p E {1,1 + 1]. It turns out that each different I requires a different construction. We 
note that there are more sophisticated constructions, when the bounds (llO.lbj) are uniform 
on compact sets, but may be not uniform in T>. The construction of |Whi36] allows the 
case of infinite derivatives. It also shows that the extended function is real analytic in the 
complement ofT>. A survey of recent developments in the problem appears in [EH]. 


Proposition 32. For given e and under the hypotheses of Theorem \TE, we have that the 
functions that produce the expansions of the normalized solution of fll.3p given by Proposi¬ 
tion^^ say Kn'^ and pin\ are continuous on Qro o.nd analytic in the interior of when 
we give the topology of Ap x C to Kn^ and fin\ 

We also have that Kn'^ and jXn'^ are C°° in the sense of Whitney in Qr^, when we give to 
the ranges the topology of Ap' x for any 0 < p' < p. 


10.5.1. Proof of Proposition Given one normalized solution, we have proved in Proposi¬ 
tion 121] that we can hnd series expansions of the normalized solutions around any point in 
Qro- These will be the Whitney derivatives. 

We note that we have produced the derivatives as polynomials in the complex variable 
£. Of course, is a multilinear function of degree j of the real variables which are the 
components of the complex number e. We will use the notation of complex powers to keep 
the notation we have used in the previous sections. 

We will now verify that the components of the expansions Kn\pin'^ verify all the properties 
in the Dehnition [301 

To verify that the functions Kn'^ , pn^ are continuous, we just observe that the Kn'^ are 
bounded in Ap/. Hence, by Montel’s theorem [SZ65j they he in a precompact set in Ap" for 
any 0 < p" < p'. When we have a function from a compact set to a compact set, continuity 
is equivalent to having a closed graph. We can check that if we have a sequence Sj —)■ e and 
Kn^'^ —t Kn\ —>■ Pn\ we Can see that the limiting functions Kn^ , pn^ satisfy (110.21) 
and it is also normalized. Using that the solution of fll.311 is unique whenever it satishes the 
normalization condition, we obtain that the solution is indeed the Kn\ 

An alternative proof of the continuity of the expansions is obtained by induction. We have 
shown that K is continuous in e and that the reducibility matrix is continuous in e in a slightly 
smaller domain in 9. The expansion is computed by recursion. It just requires to evaluate a 
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polynomial in the previous terms, multiplying by the automatic reducibility matrix, solving 
the cohomology equation, and applying other operations. Hence, by recursion, it is easy to 
see that we obtain that the terms of the expansions are continuous, provided that at every 
step we decrease (by an arbitrarily small amount) the domain of analyticity in 6 . 

We also note that Corollary |25] shows that the derivatives satisfy the Taylor estimates. 
This is the condition (llO.lbh for /c = 0. Of course, in the interior points, we obtain that 
this is indeed the derivative and, since it is a complex derivative, the function is analytic in 
a neighborhood (this gives an independent proof of the analyticity of the functions in the 
interior of Qro)- 

The only thing that remains to verify in the Dehnition [30] of Whitney derivatives is the 
condition (110.161) for k > 0. By the previous remark, we could restrict ourselves to the case 
that one of the points is close to the boundary; nevertheless, we will not need that. 

As we mentioned before, the condition for fc > 1 in Dehnition [30l unfortunately, is not 
automatic. In [dlLOOOl ILdlLlO] one can hnd examples where this fails. 

The set we are considering, however, is much simpler than many other closed sets. It is a 
connected set and indeed it has a big interior. The key property is that the sets we consider 
have the compensation property, namely for any pair of points in the domain, one can hnd a 
path joining them, whose length is less or equal than a constant times the distance between 
the points. The set Qr^ satishes a very strong form of the compensation property as shown 
in Proposition [33] below. 

Proposition 33. Given two points ei,e 2 G Gro, there is a differentiable path 7 starting in 

El, finishing in 82 and contained in Gro, indeed, contained in the interior of Gro except for 

the endpoints. Moreover I7I, the length of satisfies the ineguality 

I 7 I < 7r|ei - £ 2 ! • (10.17) 

Proof. Of course, if the domain Gro ^ ^ ^^7 convex set - the result would be 

trivial, since we could join two points £ 1,82 by a straight line which stays in the disk. 

The only thing that we have to consider is the case that some of the excluded circles from 
the boundary interrupts the straight line. 

We consider the ball B centered at the boundary of {e G C : |A(£)| <1} that passes 
through £ 1 , 82 . This ball can be easily constructed, since the center is the mediatrix between 
the points £ 1,82 and the boundary of the disk. This point is unique if we are considering A 
sufficiently small, because of the implicit function theorem. 

Now we observe that any ball centered in the boundary which intersects the segment £i —£2 
and which does not include any of the points £ 1 , £2 has to be contained in B. 

Then, to connect the two points we can move along the boundary of the disk B. Of course, 
being an arc of disk the path is bounded by tt times the distance of the points (the bound 
is saturated when the points are in opposite sides of a diameter). □ 

Using Proposition [SS] the result of Proposition [32] follows by a standard argument (see 
|LdlL10] ). Given £ 1,82 G Gro^ let 7 be a differentiable path contained in Gro^ starting in £i 
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and ending in 82 - Using the fundamental theorem of calculus, we write, as in the Lagrange 
version of Taylor’s theorem: 

_ ^{ei) ^ ^ y 

-‘'■n -‘'■n / ^ ^ ^ ^ 

i=i 

fsi fSN-i 

+ J dsi J dS 2 ---J dSN -J^p;^Kly^"'^''^'"'’^^hi{siS 2 )---{siS 2 ---SN-l) 

7'(si)7'(siS 2) • • • 7'(siS2 ■■■Sn) ■ 

Since the path 7 is contained in the interior of the domain, we have already shown that 
the terms of the Lindstedt expansions are the complex derivatives of the function and 
therefore that is analytic in the interior. Hence, it follows that the higher order derivatives 
are the derivatives of the lower order derivatives and the functions Kn^ are analytic in e. 
Hence, we obtain 

^ y 

ri -‘'■n j\ Apj ^ ^ 2 

i=i 

which is the condition (llO.lhh for Whitney differentiability. 

The proof for /i goes along the same lines and is even more elementary. □ 

Appendix A. The case of flows 

In this Section, we extend the results of the paper to the case of flows. We recall that a 
family of vector helds is conformally symplectic, if there exists a function x = x{^) with 
x(0) = 0, such that 

= X{£)^ , 

(see fll.2l) L The invariance equation for is the differential equation: 

d^Ky9)=Y^^,oKy9) , (A.l) 

where, for short, we denote the partial derivative as di_j = oj ■ de, the unknowns are K^, He- 
We look for Lagrangian tori satisfying the relation 

DKy9Y JoKy9) DKy9) = Q. 

The dehnitions [3] and 0] on the Diophantine property of the frequency need to be modihed 
as follows. 

Definition 34. Let u G r G M+. lUe define the quantity vioj] r) as 

v{uj]t) = sup \uj ■ k\~^\kY . 
fcezrf\{o} 

In the supremum above we allow infinity and we set //(a;; r) = 00 , whenever \u ■ k\ = 0. 

We say that u EMfi is Diophantine of class r and constant if 

z/(a;; r) < 00 

and we denote by Vdiv^r) the set of Diophantine vectors in M'’* of class r and constant v. 


< C\e2 — 
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We also have the following dehnition. 

Definition 35. Let x G C, ca G r G M+. We define the quantity z/(x;a;,r) as 

sup \iu ■ k + x\~^\k\~'" . 
kezfiio} 

We remark that, setting x = Xr + ^Xii then we have the following ineqnality: 

(Ku; -k + xl l/^r) >\iuJ-k + x\‘^>xl + {^-k + Xif > xl > 

from which we conclude that for Xr 7 ^ 0: 

iyiX',(^,T) < \Xr\~^ < 00 

for any r. Therefore, the critical case becomes Xr = 0, for which it could happen that 
u{x'-,uJ, t) = 00 (for example when Xi = ' k). 

We assume that the function y = x{^) is analytic in a neighborhood of e = 0 and that 
x(0) = 0; this leads us to require that y satisfies (compare with assumption HA for the 
mapping case): 

Hx x(^) = + 

for some a G Z_|_, a G C. 

In analogy to fl2.5l) and fl2.7p . we define the set Q = G{A; u, r, N, x) and its preimage under 
X, say A = A (A; u),t,N), as follows: 

g{A-,uj,T,N,x) = {eeC: iy{x{e)-,uj,T)\x{e)f^^<A} 

A{A;uj,t,N) = {xeC: i^(x;a;, r)|x(£)|^+^ < A} , (A.2) 

where A > 0, G Z_|_. Moreover, let 

gro{A-,uj,T,N,x) = g r\{e eC : l^l < ro} . (A.3) 

The equivalent of the cohomological equation fl2.8l) is given by 

d.m + A m = x{0) , (A.4) 

where X = a(^) is a known function with zero average. The solution of flA.4D can be found 
through a result extending Lemma [8], under the assumption that ca is a Diophantine vector. 

The automatic reducibility is ensured by the following result, which has been proven in 
|CCdlL13c] . 

Proposition 36. (see jCCdlL13cj ) Let N{e) = {DK{eyDK{e))-^ and let 

M{e) = [DK{e) I jDK{e)N{e)]. 

Setting A = o K, one has: 

- A(e)M(e) = m(«) 

with 

s{e) = N{e)DK{eyj(^A{e) + A{ey^DK{e)N{e). 


(A.5) 
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We now assume to start with an approximate solution {Ka,^a), satisfying the invariance 
equation flA.lj) : 

d^Ka{e)-j^^^^,oKa{e) = E{e) 

with an error term E = E{9). Then, we can prove the equivalent of Theorem [141 provided 
that the non-degeneracy condition is replaced by 


det 


7 ^ 0 , 


A S{B^)^+ A, 

yld A 2 

where S is dehned in flA.Sp . A is given by 

A = [I 1 II 2 ] = o K) 

and {B}jY solves the differential equation: 

d.{B,f + x{B,f = -{A^f . 


For completeness, we conclude by stating the main theorem for the case of vector helds. 


Theorem 37. Let AT = T'^ x B with C an open, simply connected domain with smooth 
boundary as in Theorem HH Let ca G satisfy Definition and let he a family of 
conformally symplectic vector fields with conformal factor x = x(^) Hy. We assume 

that for £ = 0, the vector field is symplectic. 

Assume that for /i = /io the flow J>o,o admits a Lagrangian invariant torus with embedding 
Kq G Ap, p > 0, satisfying 

dJ<o{0)-Ep„eoKo{e) = 0 . 


Define the following guantities: 


No{e) = 

(^DKo{9yDKo{9)) ' 

Mo{e) = 

[DKo{9)\JDKo{9)No{9)] 

Aq = 

[^10 ^^ 20 ] = —Mq ^ (VpE 0 Kq) 

A,{9) = 

0 Kq 

So{9) = 

Nq{9)DKq{9)'^ j[Aq{9) + Aq{9)^^DKq{9)Nq{9) 

Assume that the torus Kq 

satisfies the following assumption. 


HND’ Let the following non-degeneracy condition he satisfied: 

detj^^o S'o(5fe^+Aio ^ ^ 

\ 0 A 20 / 

where {Bbofi is the solution of the cohomology eguation 

d^{Bbof + x{Bbof = -{A2of . 


Then, we have the following result. 
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A) We can find a formal power series expansion Ke\ pe satisfying (lA.lj) in the sense 
that, expanding in power series and taking a truncation to order G N 


as 


k[<n] _ 




0 


[<^^1 _ 


= Ihen we have: 


O A'K«I||,, < C„|er+‘ 

for some 0 < p' < p and for some constant Cn > 0. 

B) For any N ^ fi, we can find a set Qro in (lA.3p with ro sufficiently small and for 
any 0 < p' < p we can find ■ Gro ~^p', W ■ Gro which are analytic in Gro o,nd 

extend continuously to the boundary of Gro; such that for every e G Gro satisfy exactly 
the invariance eguation: 

dujKe - oK^ = 0 . 


Moreover, the exact solution has the formal power series of part A) as an asymptotic expan¬ 
sion, in the sense that 


||A-pi-A- 

- A 


slip' 


< CN\sf^^ , 

< CN\e\^+^ . 


The proof of Theorem [37| is similar to that of Theorem fT^ and is left to the reader; it relies 
on Theorem 20 in |CCdlL13^ . 
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